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A system  for  remote  measurements  of  air-sea  fluxes 
of  momentum,  heat  and  moisture  during  moderate 
to  strong  winds 

I.  INTRODUCTION 

This  report  is  a description  of  an  internally  recording  system  for 
obtaining  estimates  of  the  fluxes  of  momentum,  sensible  heat  and  moisture 
by  the  Reynolds  flux  and  dissipation  methods.  While  these  techniques  are 
not  new  an  automatic  system  using  them  is  novel  and  so  is  the  way  the  Gill 
anemometers  are  used.  Some  of  the  techniques  should  also  be  useful  for 
acquisition  and  recording  of  other  kinds  of  geophysical  data.  Before 
describing  the  system,  a brief  review  of  the  reasons  why  such  a system 
is  required  and  of  the  existing  information  on  which  the  design  is  based 
is  qiven. 

Over  the  past  ten  years  or  so  there  have  been  many  reports  of  measurements 
of  the  momentum  flux,  some  of  sensible  heat  flux  and  a few  of  moisture 
flux.  While  there  is  considerable  scatter  in  each  data  set  for  a variety 
of  possible  reasons,  a reasonably  consistent  picture  for  the  momentum 
flux  seems  to  have  emerged  (Hidey,  1972;  Phillips,  1972;  Pond,  1975; 
Stewart,  1974)  over  the  range  of  wind  speeds  for  which  good  measurements 
have  been  obtained  (up  to  about  15  m/s).  It  appears  that  this  flux  can 
be  parameterized  in  terms  of  the  wind  speed  squared  with  a non-dimensional 
coefficient  - the  drag  coefficient,  C^.  For  the  other  fluxes  there  are 
not  so  many  observations  but  it  appears  that  similar  parameterizations  may 
be  possible  for  them  as  well  (Pond,  et  al.,  1974).  Even  for  the  momentum 
flux  there  are  insufficient  data  to  try  to  pick  out  particular  sources  of 
variation  within  the  real  statistical  scatter  such  as  those  caused  by 
stability,  stage  of  development  of  the  wave  field,  presence  of  a frontal 
system,  and  lack  of  complete  stationarity.  There  are  also  two  schools  of 
thought  on  whether  or  not  the  drag  coefficient  depends  on  wind  speed. 
Phillips  (1972)  and  Stewart  (1974)  suggest  that,  in  view  of  the  scatter,  a 
constant  value  is  reasonable  over  the  existing  range  of  measurements  while 
Smith  and  Banke  (1975)  and  Garrett  (1977)  among  others  suggest  some 
increase  with  wind  speed.  Because  of  the  small  amount  of  data,  such  an 
increase  could  be  found  in  some  sets  of  data  or  could  be  missed  in  other 
sets  due  to  other  factors.  For  example,  Denman  and  Miyake  (1973)  suggest 
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that  during  conditions  of  increasing  wind  velocity  the  drag  coefficient 
is  higher  because  the  wave  field  is  not  in  equilibrium  with  the  wind.  In 
decreasing  winds  one  might  expect  the  opposite  effect.  In  any  case,  the 
trend  is  small  enough  that  it  probably  would  have  little  effect  on  computed 
average  stress  values,  provided  that  winds  over  15  m/s  do  not  make  important 
contributions,  although  calculations  of  the  curl  of  the  stress  might  be 
affected  by  the  difference.  If  the  present  data  gave  results  to  20  m/s  it 
would  be  possible  to  make  more  definite  statements.  For  the  stress  vector 
averaged  over  a month  or  so  there  are  few,  if  any,  places  on  the  world  ocean 
where  winds  greater  than  20  m/s  make  important  contributions  because  they 
occur  too  infrequently. 

Clearly  there  are  a number  of  reasons  why  one  would  like  to  collect  more 
data  to  calculate  values  of  Cg.  In  addition,  if  one  is  interested  in 
following  the  time  history  of  the  stress  and  its  effects,  then  it  becomes 
very  important  to  extend  the  measurements  to  the  higher  wind  speed  range. 

Such  measurements  should  settle  the  controversy  over  the  trend  and  give 
a much  better  measure  of  it.  More  measurements  of  the  fluxes  of  sensible 
heat  and  moisture  and  of  their  parameterization  coefficients  are  needed  at 
all  wind  speeds,  particularly  high  ones. 

Such  measurements  have  not  been  obtained  so  far  because  of  the  experimental 
difficulties.  However,  it  should  be  possible,  because  of  the  development 
of  techniques  for  obtaining  the  necessary  data,  because  of  the  knowledge 
of  the  nature  of  quantities  now  gained,  and  because  of  developments  in 
modern  electronics,  to  design  a system  which  will  allow  measurements  at 
higher  speeds.  One  of  the  serious  difficulties  is  to  find  a suitable  plat- 
form from  which  such  measurements  can  be  made  under  open  sea  conditions. 

For  the  'best'  technique  (the  Reynolds  flux  method)  a stable  platform  is 
required.  In  addition,  one  would  like  to  develop  methods  which  would  work 
on  more  readily  available  platforms  such  as  weatherships  and  other  ships 
of  opportunity  (e.g.  the  dissipation  method,  provided  comparison  with  the 
Reynolds  flux  method  shows  that  it  is  satisfactory).  Internal  recording 
would  seem  to  be  desirable  because  any  suitable  stable  platforms  would 

probably  be  unmanned  (for  reasons  of  safety,  among  others)  and  because  4 

for  a ship  of  opportunity  program  it  is  much  more  efficient  not  to  have 
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to  send  a team  of  observers.  Some  compromises  would  have  to  be  made  with 
this  choice.  Telemetry  is  an  al ternati ve, but  would  probably  not  be  suitable 
for  a ship  of  opportunity  program,  increases  the  power  requirements  sub- 
stantially, is  more  complicated  and  has  range  limitations. 

Finally,  it  is  worth  noting  that  to  be  useful  for  other  purposes,  such  as 
calculating  fluxes  over  the  world  ocean  or  parts  thereof  for  modeling  and 
predicting  the  behaviour  of  the  ocean  and  atmosphere,  more  or  less  direct 
flux  observations  must  be  used  to  obtain  parameter izat ions  that  estimate  the 
fluxes  from  more  easily  measured  and  readily  available  variables.  More  or 
less  direct  flux  observations  are  simply  too  difficult  and  too  costly  to 
be  made  over  wide  areas  on  a routine  basis.  Generally,  the  direct  measure- 
ments have  been  used  to  examine  and  evaluate  coefficients,  such  as  Cp.  in 
the  bulk  aerodynamic  parameteri zat ion  (given  in  the  next  section).  Even 
the  data  required  to  estimate  fluxes  using  this  parameterization  are  not 
as  available  as  one  would  like.  Thus  another  part  of  air-sea  interaction 
research  to  be  pursued  is  a search  for  parameteri zatlons  in  terms  of  readily 
available  data,  e.g.  wind  roses  and  long-term  average  sea  surface  and  air 
temperatures  and  humidities  for  long-term  average  fluxes  or  surface  pressure 
naps  and  temperature  and  humidity  values  to  look  at  time  dependent  processes 
(Fissel,  Miyake  and  Pond,  1977). 

1 1 • THEORY  AND  EXISTING  RESULTS 

Only  an  outline  will  be  given  here  as  more  detailed  treatments  are  available 
in  many  sources  (e.g.  Roll,  1965;  Lumley  and  Panofsky,  1964;  Monin  and 
Yaglom,  1965,  1967;  Kraus,  1972;  Burling  and  Stewart,  1967;  Busch  and 
Panofsky,  1963;  Miyake,  Stewart  and  Burling,  1970;  Pond,  1975;  Busch,  1977). 

a)  Reynolds  flux  method  - 

This  is  the  most  direct  method  of  estimating  the  fluxes  and,  with  the 
development  of  techniques  to  make  the  necessary  measurements,  is  generally 
taken  to  be  the  standard  to  which  other  techniques  are  compared  directly, 
or  Indirectly  through  comparison  of  parameterization  coefficients.  It 
depends  on  measuring  the  covariance  between  w,  the  vertical  component 
of  the  velocity,  and  U,  T and  q,  the  downstream  velocity  component, 
temperature  and  absolute  humidity,  respectively.  Since  the  average  of 
w(w)  must  be  zero  (Burling  and  Stewart,  1967),  the  total  values  (U,T,q) 


f 


or  the  fluctuations  (u  = U - U , T'  * T - T,  q‘  = q - q)  may  be  used . Thus  we 
have: 

Stress  magnitude,  t * -p  Dw  » -p  uw  and  acts  in  the  direction 

of  the  wind 

Sensible  heat  flux,  = pC  wT  = pC  wT1 

^ r P 

Moisture  flux,  E = wq  = wq'  = (latent  heat  flux,  H^VL 

where  p is  the  air  density  (whose  fluctuations  are  small  enough  to  be 
neglected  in  the  calculations),  is  the  specific  heat  of  air  at  constant 
pressure,  and  L is  the  latent  heat.  In  practice,  one  must  nake  sure  that 
w does  not  have  a mean  or  large  trends  introduced  spuriously  by  the  instru- 
ments. Generally,  one  should  use  the  variations  about  the  mean  of  the  record 
used.  Whether  to  detrend  or  not  is  not  clear-cut.  Personally,  unless  we 
suspect  considerable  instrument  drift,  we  prefer  not  to  detrend.  Other 
groups  always  do.  Either  method  may  introduce  some  unmeasurable  but 
hopefully  small  errors  in  the  results. 


It  is  useful  to  examine  the  cospectra  of  these  quantities  to  ensure  that 
all  contributions  are  being  included.  Contributions  occur  over  the  normal- 
ized frequency  range,  0.001  < n v 10,  where  n ■ fZ/U,  f is  the  frequency 
(in  Hz),  Z is  the  observation  height  (in  m)  and  1)  is  the  magnitude  of  the 
mean  velocity  (in  m/sec).  Thus  one  should  work  with  a record  length 
T 1000Z/U.  In  practice,  a few  such  records  should  be  averaged  together 
to  get  reasonable  statistical  stability  of  the  low  frequency  estimates. 

E.g.,  for  Z/U  * 1,  four  15-minute  records  or  an  hour  of  data  should  be 
adequate.  Because  of  the  large  frequency  ranoe  it  is  usual  to  plot  log  f 
or  loq  n.  To  obtain  a plot  whose  integral  is  proportional  to  the  covariance 
and  such  that  the  relative  contributions  from  different  frequency  bands 
are  clearly  shown,  one  then  plots  frequency  times  cospectrum,  C0,  since  then 
2.3f  C0  d ( 1 og  f)  - Co  df,  and  the  J"  a covariance.  In  the  following, 
when  we  refer  to  spectra  and  cospectra  we  shall  mean  frequency  x these 
variables  unless  otherwise  noted.  These  cospectra  normalized  according  to 
the  Monin-Oboukhov  similarity  theory  and  plotted  aoainst  log  n should  be 
functions  of  stability.  However,  for  the  stabilities  usually  encountered 
at  sea  and  within  the  natural  scatter,  stability  effects  can  probably  be 
ignored.  The  peaks  of  the  cospectra  occur  for  0.01  ^ n s 0.1  in  most  cases 
although  sometimes  the  peak  of  the  wT  cospectrum  occurs  at  values  of  n 
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up  to  about  0.5.  The  contribution  to  the  cospectra  from  n > 1 is  only  a 
few  percent  of  the  total  except  for  these  unusual  cases  of  the  wT  cospectrum 
where  the  contribution  may  be  up  to  20*  of  the  total.  It  is  now  thought  that 
these  cases  of  wT  cospectra  with  large  contributions  at  larger  n are 
spurious  results  caused  by  moisture  films  on  the  temperature  sensors  which 
introduce  large  errors  when  the  latent  heat  flux  is  large  compared  to  the 
sensible  heat  flux  (Friehe,  et  al.,  1976). 

Measurements  of  this  type  are  made  some  distance  above  the  surface  but  it 
is  not  difficult  to  show  that  the  fluxes  are  nearly  constant  for  the  first 
few  tens  of  metres  (Lumley  and  Panofsky,  1964).  The  thickness  of  this 
'constant  flux'  layer  generally  should  increase  with  the  wind  speed. 

Assuming  that  the  normalized  cospectra  will  not  be  grossly  different  at 
high  wind  speeds,  this  information  allows  one  to  work  out  a suitable 
sampling  scheme  in  terms  of  rate  and  record  length. 

While  it  is  perhaps  the  best  method,  the  Reynolds  flux  technique  does  have 
some  limitations.  In  order  to  obtain  the  cospectra  to  ensure  that  the  co- 
variances  are  correct  one  must  sample  at  a fairly  fast  rate  for  a fairly 
long  time,  so  a rather  large  amount  of  data  must  be  recorded  to  get  a single 
Reynolds  flux  estimate.  The  method  requires  a known,  preferably  fixed, 
instrument  orientation.  In  principle,  it  is  possible  to  use  a platform 
such  as  a surface  ship,  measure  the  motion  and  make  corrections  point  by 
point.  Such  an  approach  makes  the  measurement  program  much  more  compli- 
cated and  thus  makes  it  even  more  difficult  to  obtain  measurements  one  can 
be  confident  in.  As  the  problem  at  high  speeds  is  difficult  enough  on  a 
stable  platform  it  seems  best  to  use  such  a platform,  at  least  initially. 

One  might  also  look  for  and  test  other  methods  at  the  same  time  which 
could  be  used  on  mobile  platforms  and  have  reduced  recording  requirements, 
b)  Pulk  aerodynamic  parameter) zat ions 


In  this  method  one  assumes  that  the  fluxes  may  be  parameterized  in  the 
followinci  wav: 


'•■•V-W 


1 


6 

|t I - pCDU2 

HS  * pCp  CT  °*T 

Hl  3 Cq  UAq  • L 

where  U is  the  magnitude  of  the  vector  mean  wind,  aT  = sea  surface  - air 
temperature,  and  Aq  = sea  surface-air  humidity  and  0,  air  temperature  and 
humidity  are  measured  at  or  reduced  to  some  reference  height,  commonly  10  m; 
Cq,  Ct  and  Cq  are  non-dimensional  coefficients.  Cq  is  called  the  drag  co- 
efficient; Cy  and  Cq  are  sometimes  referred  to  as  the  Stanton  and  Dalton 
numbers,  respectively.  U,  AT  and  Aq  should  be  averages  over  times  of 
order  an  hour.  Generally,  sea  surface  temperatures  (and  hence  humidity  at 
the  surface)  are  based  on  occasional  spot  readings,  usually  'bucket'  samples. 
Spot  readings  or  fairly  short  averages  (the  usual  ship  observations)  may  be 
used  with  some  additional  statistical  scatter  since  the  means  of  U,  aT  and 
Aq  are  generally  much  larger  than  the  root  mean  square  (rms)  fluctuations. 

To  calculate  the  parameterization  coefficients,  ||J|,  the  magnitude  of 
vector  mean  wind,  or  the  mean  speed  over  the  record  length  of  about  an 
hour,  is  used  depending  on  the  data  set.  Normally,  for  periods  of  an 
hour  or  so  the  mean  speed  is  less  than  1%  higher  than  the  magnitude  of 
the  vector  mean  (see  Appendix)  and  the  differences  can  be  ignored.  One 
can  examine  Cq,  Cy,  and  Cq  to  see  whether  they  are  reasonably  constant  or 
have  measurable  dependence  on  wind  speed,  stability,  wave  development,  etc. 
One  can  also  use  these  formulae  to  estimate  the  ratios  of  scaling  amplitudes 
to  mean  values.  For  this  purpose,  to  set  design  criteria  without  making 
any  claim  to  exactness,  we  take  CQ  'v  Cy  ^ Cq  ^ ( 1/25) 2 = 1.6 x 10'3.  We  put 

C0n2  ■ “.2 

Cy  UaT  s X U*  T* 

Cq  UAq  = k u*  q* 

L 

u*  = | uw | is  called  the  friction  velocity  and  is  the  velocity  scale  of  the 
Monin-Oboukhov  similarity  theory,  k is  von  Karman's  constant  (*0.4)  and  is 
included  here  in  the  definition  of  the  scaling  parameters  T*  and  q*  as  is 
usually  (but  not  always)  done.  Then 

u*  * 1/25  IT,  T*  s 1/10  aT  and  q*  a 1/10  Aq. 
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c)  Statistical  quantities 

Here  we  want  to  relate  amplitudes  of  the  fluctuations  to  the  scales 

u*,  T*  and  q*  and  hence  to  the  mean  values,  II,  aT  and  Aq.  Again,  we 

present  values  that  do  not  necessarily  provide  the  final  answer  but  are 

useful  for  design,  e.g.  for  picking  amplifier  gains.  We  tend  to  pick 

numbers  near  the  high  end  of  existing  data  to  avoid  overload  problems. 

We  look  at  rms  values  oa  = (a^)'2  where  a nay  be  u,  w,  T'  or  q'. 

aw/u*  1.2 -1.5  so  we  take  aw  =0.06  U,  perhaps  a slightly  high  estimate. 

a,,/u*  ^ 2.5  so  we  use  a =0.1  U.  o /u*  values  show  considerable  scatter 
u u u 

which  is  associated  with  the  fact  that  low  frequencies  dominate  the  contri- 
butions tolls'.  Thus,  the  record  length  and  detrending  can  change  this 
value  considerably.  These  values  should  be  representative  if  the  low 
cut  off  normalized  frequency  n -v  0.001  to  0.002.  The  correlation 

between  u and  w,  -u £/a  a for  a = 0.06  U,  a = 0.1  U,  and  u*  = 1/25  U 

u w w u 

is  -0.27  which  is  close  to  observed  values.  All  these  ratios  may  have  some 
dependence  on  stability  but  it  is  too  small  to  be  well  established  experi- 
mentally and  can  be  ignored  for  design  purposes  over  the  expected  range. 

Oj/ T*  and  a^/q*  seem  to  have  stronger  stability  dependence  although  it  is 
not  well  established  either.  Oy/T * ^ oq/q*  % % to  1 so  we  take  ay  -v  1/10  aT 

and  0^  ^ 1/10  Aq  which  we  expect  to  be  reasonable  upper  limits. 

Clearly,  since  none  of  these  numbers  are  very  exact  some  leeway  must  be 
allowed  in  our  design.  In  addition,  if  CD  does  increase  with  wind  speed 
some  other  things  must  increase  too.  Over  land  C^'s  are  usually  consider- 
ably larger  and  fairly  independent  of  IT  provided  the  geometry  of  the  surface 
is  independent  of  U,  but  a.,/u*  and  o ,/u*  and  the  uw  correlation  remain  about 
the  same.  If  Cq  is  higher  over  water  at  higher  wind  speeds,  it  may  be  simply 
because  of  a higher  apparent  surface  roughness  with  all  the  properties 
scaled  according  to  the  Monin-Oboukhov  theory  remaining  the  same  as  they 
do  for  existing  measurements  over  both  land  and  water.  One  must  allow 

for  such  a possibility  and  higher  values  of  a /U  and  a / U in  any  design. 

u w 

d)  The  dissipation  method 

This  method  of  estimating  the  fluxes  does  not  require  the  vertical 
component  of  the  velocity,  w,  and  is  thus  more  suited  for  use  on  a ship 
or  mobile  platform  than  is  the  eday  flux  method.  It  is  based  on  the  idea 
that  the  production  of  turbulent  fluctuations  in  the  velocity  or  scalar 
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fields  is  balanced  locally  by  molecular  dissipation.  In  unstable  conditions 
this  balance  in  the  kinetic  energy  equation  includes  a buoyant  production 
term.  This  unstable  case  and  possible  approaches  by  which  one  may  be  able 
to  allow  for  it  are  discussed  by  Pond  et  al.  (1971)  or  Rusch  (1977). 

Here,  for  brevity,  we  shall  discuss  the  neutral  case  only.  Over  the  ocean 
at  moderate  to  strong  winds  the  stability  is  almost  always  near  neutral. 
Exceptions  nay  occur  over  regions  such  as  the  Gulf  Stream  or  Kuroshio  with 
cold  air  outbreaks  from  the  continents.  There  are  other  terms  in  the  energy 
equations  but  the  available  data  show  that  producti<*s  dissipation  for 
near  neutral  to  moderately  unstable  cases.  Ue  expect  this  balance  to 
continue  to  hold  at  high  wind  speeds  (because  according  to  the  Monin- 
Oboukhov  similarity  theory,  which  rationalizes  existing  data  over  both 
sea  and  land,  such  relations  are  not  wind  speed  dependent)  but  an  essential 
part  of  our  measurement  program  is  to  verify  this  expectation.  If  we  find 
the  balance  does  not  hold,  then  we  shall  attempt  to  find  an  empirical 
correction  which  allows  us  to  estimate  the  fluxes.  (Our  preliminary  analysis 
of  some  of  our  data  shows  that  obtained  from  the  dissipation  method  is 
statistically  similar  to  Cp  obtained  by  the  Reynolds  flux  method.)  For 
near  neutral  stability  and  stationary  and  horizontally  homogeneous  flow, 
the  balance  is: 

5 l 

dissipation 

« N'* 

where  > and  >’  are  used  to  represent  the  mean  and  fluctuating  part, 
respectively,  of  any  scalar  field. 

The  mean  gradients  are  given  by  the  Monin-Oboukhov  similarity  theory 
dU  , u* 

37  7[ 

af  ■ ¥ 

Thus,  -uw  * u**  = (kiZ)"  3 

|wT’  I • |mi*uI  a UZ)2/3  nX1/6. 
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The  sign  of  w>'  is  given  by  the  sign  of  the  air-sea  surface  difference  in 
y.  These  equations  allow  estimates  of  the  fluxes  to  be  made  from  measure- 
ments of  e and  Ny. 


Assuming  that  Kolmogoroff  theory  is  applicable,  the  one  dimensional  spectra 
of  downstream  velocity  (4^  and  scalar  fluctuations  (4^)  for  values  of  k , 
the  downstream  component  of  the  radian  wave  number,  in  the  inertial  subrange 


are  given  by: 


K'  c2'3  k^3 
B ' N e'1/3  k"5/3 

V v 


K’  and  ’ are  Kolmogoroff  constants  (0.55  and  0.8,  Paquin  and  Pond,  1971), 
and  the  observed  frequency  f and  mean  wind  U give  k using  Taylor's 


hypothesis  (k  = 2tt f / U) . 


The  dissipation  of  mechanical  energy,  e , is  thus  obtained  solely  from 
measurements  of  the  downstream  velocity  component,  U,  provided  that  fre- 
quencies extending  into  the  -5/3  range  of  <j>  are  observed.  In  practice, 
this  measurement  is  not  too  difficult,  as  a normalized  frequency  of  n > 0.2 
should  be  adequate.  The  -5/3  range,  if  isotropy  is  really  required,  should 
not  start  until  n > 1,  but  empirically  it  extends  back  to  n ^ 0.2  or  even 
0.1  for  the  downstream  velocity  and  for  temperature  and  humidity,  which  is 
fortunate  because  relatively  simple  robust  sensors  which  respond  fully  past 
n = 1 are  not  available.  Considerable  platform  motion  is  acceptable. 

TSupposed  the  measured  component  is  U cos  4 + w sin  <5  + UR,  where  UR  is  the  addi- 
tional velocity  caused  by  the  fluctuating  platform  motion  and  6 is  the  tilt  into 
the  wind  (see  Appendix).  Provided  we  make  our  measurements  outside  the  fre- 
quency range  of  the  platform  motion,  UR  can  be  ignored.  Then  the  measured 

spectrum  is  4 cos^T  + 4 sTnTs  + 4W  sin-  5 , assuming  >5  is  not  correlated  with 

u uw  " 

u and  w which  is  reasonable.  Mow  |6|  should  be  fairly  small  so  we  get  approximately 
4(1-  7^)  + 2? 4 + 52  4 . With  \7\  <*  5°  and  0.  = 10°,  a fairly  extreme  case, 

1 6 1 < 0.1,62  <_  0.05.  Further,  4=4  and  |4  J < 4u  , so  there  may  be  a 
random  error  of  up  to  10?  in  this  extreme  case.  Almost  all  of  this  error  is 
associated  with  7 which  one  would  hope  would  not  be  as  large  as  5°  in  r’ze. 

|4  |/4  depends  quite  strongly  on  n;  for  n > 0.3  where  one  would  try  to  work 

this  ratio  < and  the  error  < 5?.l  Measurement  of  a scalar  variable  in  the 
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-5/3  range  of  its  spectrum,  gives  the  dissipation  of  scalar  fluctua- 
tions Ny . In  principle,  only  a single  discrete  value  from  each  of  the 
spectra,  averaged  over  a few  minutes,  is  required.  The  recording  require- 
ments are  greatly  reduced  compared  to  the  Reynolds  flux  method.  Time 
histories  of  the  fluxes  may  be  compiled.  In  practice,  we  use  three  fre- 
quency bands  for  each  variable  as  will  be  discussed  later. 

The  dissipation  method  has  been  shown  to  give  values  of  the  momentum  flux 
comparable  to  those  of  the  eddy  flux  method  in  several  studies  (e.g.  Pond 
et  al.,  1971;  Smith  and  Banke,  1975),  mainly  in  light  to  moderate  winds. 

Only  a few  comparisons  of  moisture  and  sensible  heat  fluxes  from  both 
Reynolds  flux  and  dissipation  methods  have  been  made  but  the  method 
appears  feasible  (Pond  et  al.,  1971,  1974). 

III.  SENSORS 

For  work  at  moderate  to  high  wind  speeds  (5  ms-1  or  10  knots  to  hopefully 

30  - 40  ms-1  or  60-80  knots)  sensors  must  be  chosen  which  are  reasonably  robust. 

For  operation  on  a remote  tower,  power  is  likely  to  be  limiting  so  a low  power 
drain  is  advantageous.  For  both  tower  operations  and  ship  of  opportunity 
operations  reliable  operation  without  servicing  for  extended  periods  is 
essential.  It  is  also  vital  to  be  able  to  tell  whether  the  data  which 
are  obtained  are  reliable.  Thus  it  is  useful  to  measure  the  same  variable 
in  more  than  one  way.  Because  under  high  wind  conditions  damage  to  the 
sensors  is  inevitable,  it  is  helpful  to  use  commercially  available  components, 
preferably  of  moderate  cost. 

Some  compromises  may  have  to  be  made  in  attempting  to  satisfy  all  these 
requirements,  particularly  those  of  calibration  stability  and  robustness, 
while  still  having  sufficient  sensitivity  to  measure  the  higher  frequency 
contributions  to  the  Reynolds  fluxes.  Working  at  high  wind  speeds  has  the 
advantage  of  improving  the  response  if  it  is  convection  limited,  e.g. 
temperature  and  some  humidity  sensors.  Non-linear  response  of  mechanical 
wind  measuring  equipment  which  occurs  under  light  winds  may  also  be  ignored. 
Also,  the  nearly  constant  flux  layer  is  thicker  so  one  can  work  at  greater 
heights.  The  sensors  must  be  put  fairly  high  to  avoid  waves  and  heavy 
spray  and  distortion  of  the  air  flow  by  the  platform  (e.g.  a ship)  in  any 


case. 
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Ideally,  one  would  like  response  to  a normalized  frequency  n = fZ/U  of  10 
but  3 is  adequate.  Response  3 db  down  at  n = 1 requires  fairly  small 
corrections  (v  20%)  and  can  be  tolerated  if  the  response  is  reasonably 
well  known.  The  heiqhts  we  have  used  are  about  13  to  20  m for  open  ocean 
data  with  some  intercomparison  data  taken  at  about  10  m.  To  get  an  indi- 
cation of  the  response  required,  take  Z/U  v 1/2  s,  then  n = 3 at  f = 6 Hz. 
If  the  response  occurs  as  a distance  constant  limitation  as  in  a mechanical 
wind  sensor,  then  the  3 db  down  frequency  will  be  LP(2itD)  where  0 is  the 
distance  constant.  To  get  n = 3 requires  D = Z/6n  or  about  1/2  m for 
Z = 10  m and  1 m for  Z = 20  m. 

a)  Temperature 

To  measure  air  temperature  and  its  fluctuations  we  use  glass-coated 
microbead  thermistors  (Victory  Engineering  Corp.).  Their  response  is  3 db 
down  at  a few  Hz  at  moderate  wind  speeds  and  improves  with  wind  speed  to 
about  30  Hz  at  70  ms"1  (Miyake  et  al.,  1970)  and  should  be  sufficient. 
(Better  response  can  be  achieved  with  resistance  wire  sensors  but  they 
are  not  as  strong,  require  much  greater  ampl i fication,  which  is  more  likely 
to  cause  problems  , and  are  likely  to  be  more  subject  to  corrosion).  The 
microbead  is  operated  in  a bridge  circuit  with  a micropower  operational 
amplifier  as  detector.  By  suitable  choice  of  a resistor  in  series  with  the 
thermistor,  the  non-linearity  of  the  bridge  and  thermistor  can  be  balanced 
so  that  the  output  is  linear  within  measurement  error  over  a range  of  about 
25  C°.  The  current  through  the  bead  is  kept  small  enough  that  velocity 
sensitivity,  because  of  variation  in  self-heating  with  flow  speed,  is 
negligible  above  1-2  ms-1. 

A glass  rod  thermistor  potted  in  soft  epoxy  is  used  to  measure  sea  temper- 
ature a few  metres  below  the  surface  when  working  on  a tower.  During  our 
first  major  field  operation  which  was  on  the  beach  at  Sable  Island,  Nova 
Scotia,  the  sea  temperature  thermistor  was  put  next  to  the  microbead. 
Comparison  of  the  two  temperature  signals  indicated  erroneous  temperatures 
from  the  microbead  after  several  days.  We  believe  that  the  differences 
were  caused  by  a saline  moisture  film  across  the  bead  and  its  leads.  To 
avoid  this  problem  and  corrosion  which  eventually  causes  failure,  we  now 
insulate  the  leads  right  up  to  the  glass  coating  using  thinned-down 
transformer  varnish  (glyptol).  The  coating  near  the  bead  is  done  with 
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a single  hair  brush  under  a microscope  and  tests  have  shown  that  complete 
insulation  is  achieved.  We  also  added  a sealed  glass  rod  thermistor  next 
to  the  microbead  to  provide  a continuous  check  on  it.  With  the  coating  the 
microbeads  seem  to  be  very  reliable  until  they  get  broken,  which  happens 
very  occasionally  in  spite  of  the  protective  cover  we  put  them  in.  In  our 
present  field  work  on  board  the  CCGS  Quadra,  one  of  the  two  ocean  weather 
station  PAPA  ships,  we  have  two  microbeads  to  help  ensure  getting  temperature 
fluctuation  data. 


The  accuracy  of  our  temperature  measurements  is  about  ± 0.2  C°;  differences 
between  sensors  are  accurate  to  ± 0.1  C°  or  better.  The  gains  of  the  bridge 
amplifier  were  selected  to  give  a full-scale  range  of  20  - 25  C°  for  the  ± 5v 
input  range  of  the  recording  system.  The  midpoint  of  the  range  is  adjustable 
with  a variable  resistor  in  the  other  arm  of  the  bridge.  The  recording 
system  has  12-bit  resolution  so  for  temperature  our  resolution  is  ± 0.006  C°. 
The  amplifiers  were  also  designed  co  act  as  single-pole  low-pass  filters 
3 db  down  at  about  15  Hz,  to  reduce  noise  outside  the  frequency  range  of 
interest. 

b)  Humidity 

In  spite  of  trying  several  possibilities  we  have  not  found  a humidity 
sensor  with  adequate  frequency  response  which  is  suitable  for  unattended 
operation  for  reasonable  periods  of  time.  We  initially  rejected  the 
a-Lyman  humidiometer  (Electromagnetic  Research  Corporation)  because  of 
relatively  high  power  requirements  (the  latest  model  is  much  better  in 
this  respect) , limited  source  tube  life,  and  calibration  drift  due  to 
window  contamination  and  source  tube  aging.  The  a-Lyman  humidiometer 
does  have  good  frequency  response  (much  better  than  any  other  humidity 
sensor  we  know  of).  It  is  limited  by  the  path  length  between  the  tubes 
(about  1 cm)  or  more  likely  by  the  flushing  of  the  protective  cover  one 
puts  over  it.  Based  on  a local  field  test  over  land  at  moderate  wind 
speed  and  comparison  with  the  a-Lyman  humidiometer,  we  found  two  other 
humidity  sensors  to  try  which  appeared  to  have  just  adequate  frequency 
response  - an  aluminum  oxide  sensor  (Panametrics,  Inc.)  and  a Brady  array 
(Thunder  Scientific  Corp.).  Both  have  low  power  requirements  - a few  mA. 

In  local  field  trials  of  our  system  it  still  appeared  that  they  would  be 
possible  to  use,  but  on  Sable  Island  where  there  is  more  salt  in  the  air 
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because  of  a much  higher  (and  more  oceanic)  surface  salinity  we  found  the 
AIO2  sensor  to  have  rapid  calibration  drift  and  a limited  lifetime.  Un- 
fortunately it  is  not  suitable  for  use  in  a marine  environment.  The  Brady 
array  showed  drift  too,  but  it  was  not  as  severe.  A second  Brady  array 
with  a protective  filter  was  added  to  replace  the  AIO2  sensor  in  an  attempt 
to  provide  a low-frequency  calibration  of  the  unfiltered  Brady  array.  The 
calibration  drift  of  the  filtered  array  is  much  reduced  but  the  analysis 
of  the  Sable  Island  data  shows  that  the  response  of  the  unfiltered  array 
is  much  poorer  than  we  hoped.  The  poor  response  may  be  a salt  contamination 
problem  too.  In  any  case,  it  appears  that  the  Brady  array  is  not  suitable 
for  measuring  humidity  fluctuations  over  the  required  range  in  an  oceanic 
environment. 

For  our  measurements  on  the  CC6S  Quadra  (which  occupies  weather  station 
PAPA  half  the  time)  we  obtained  a dew  point  system  (Cambridge  Systems 
Model  2000)  to  provide  a low-frequency  calibration  of  an  a-Lyman  humidiom- 
eter.  It  also  provides  another  air  temperature  to  check  our  other  measure- 
ments. It  seems  to  work  for  a month  or  so  before  salt  contamination  makes 
cleaning  of  the  dew  point  mirror  necessary.  Unfortunately,  serious  con- 
tamination of  the  a-Lyman  windows  is  quite  rapid,  within  a very  few  days 
at  best.  We  plan  to  use  the  dew  point  system  and  a-Lyman  humidiometer 
during  JASIN  (Joint  Air-Sea  Interaction  Experiment)  in  the  North  Atlantic 
in  1978  when  we  shall  be  on  board  and  can  do  the  necessary  servicing  in 
another  attempt  to  measure  and  parameterize  the  moisture  flux. 

c)  Velocity 

There  are  a large  variety  of  sensors  to  consider  for  this  measurement. 
However,  the  conditions  of  being  able  to  operate  in  rain,  spray  and  perhaps 
snow  or  hail,  and  having  good  calibration  stability,  eliminate  most  of  the 
possible  sensors.  We  choose  to  base  our  velocity  system  on  the  GILL 
propeller-vane  system  (R.M.  Young  Co.)  based  on  the  experience  of  others 
in  our  Institute  (e.g.  Denman  and  Miyake,  1973)  among  many  others.  It  is 
a mechanical  system  which  has  good  calibration  stability  not  only  with  a 
particular  propeller  but  also  from  propeller  to  propeller,  since  the 
geometry  is  accurately  maintained  (even  rather  severe  ablation  of  the 
propellers  by  rain,  spray  and  hail  does  not  affect  the  calibration 
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measurably),  which  is  simple  and  easy  to  service,  which  has  very  low  power 
requirements  (for  the  preamplifiers  to  give  the  tacogenerators  the  desired 
output  range  and  to  supply  voltage  to  the  direction  potentiometer),  and 
which  works  reliably  for  fairly  long  periods  for  moderate  to  strong  winds. 
(For  light  winds  the  friction  limits  of  a mechanical  system  become  important 
but  they  are  not  of  concern  to  us.)  Furthermore,  the  sensitivity  of  the 
system  is  quite  good.  Quoted  values  of  the  distance  constant  are  about  a 
metre  which  is  adequate  if  the  observation  height  is  large  enough  and  if 
corrections  are  made  when  necessary.  (A  method  of  evaluating  the  response 
will  be  discussed  in  the  next  sub-section.) 

These  sensors  have  been  used  in  a variety  of  ways.  One  is  to  use  three 
propellers  with  their  axes  at  right  angles  to  give  three  components.  Another 
is  to  use  a propeller  vane  for  the  horizontal  components  and  a propeller  with 
a vertical  axis  for  the  vertical  component.  The  problem  with  this  second 
system  is  that  the  vertical  component  measurement  is  affected  by  non-linear 
response  associated  with  friction  and  the  driving  of  the  tacogenerator  for 
vertical  components  near  zero.  These  problems  occur  for  a velocity 
component  less  than  about  ± 1 ms"1,  although  some  correction  can  be  made 
by  adjusting  the  calibration  coefficient. 

To  avoid  this  problem  we  modified  the  propeller-vane  system  by  extending  the 
vertical  supporting  shaft  and  adding  a second  propeller  with  an  axis  about 
30°  from  the  vertical.  We  have  made  this  modification  to  some  of  the 
factory  systems;  we  have  also  had  the  factory  make  a long  vertical  shaft 
system  for  us  and  added  the  attachment  for  the  second  propeller.  The 
standard  system  of  slip  rings  is  already  adequate  with  four  slip  rings  with 
dual  contacts  (presumably  for  compatibility  with  a propeller  bivane  system, 
also  manufactured);  these  are  used  in  pairs  for  the  two  propellers.  By 
using  a 30°  axis  propeller  we  introduce  some  of  the  horizontal  component 
(approximately  %)  into  this  measurement.  Thus  this  propeller  always  turns 
fairly  rapidly,  eliminating  the  region  of  stall  and  non-linear  response 
for  wind  speeds  larger  than  4-5  ms-1. 

A photograph  of  the  twin  propel ler-vane  system  is  shown  in  Figure  1.  The 
protective  cover  and  radiation  shield  for  the  temperature  (and  Alf^  or 


sensor  level 


a)  Gil1  twin  propel ler- vane  anemometer  (propeller  diameter  0.19  m) 

b)  Enclosure  for  temperature  and  humidity  ( A 1 0 . and  Brady)  sensors 

c)  Cross-section  of  flov  path  throuqh  the  enclosure. 
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Brady  humidity  sensors  when  they  were  used)  is  also  shown  in  Figure  1. 

The  cylinder  on  top  houses  the  temperature  bridges  and  humidity  preamplifier 
circuits.  The  anemometer  preamplifier  circuits  are  inside  the  housing  at 
bottom  of  the  anemometer.  The  preamplifiers  have  full-scale  ranges  (±  5v) 
of  ± 180°  for  direction,  of  50  ms-1  for  the  horizontal  axis  propeller  and 
35  ms-1  for  the  tilted  propeller;  they  also  act  as  single-pole  filters  3 db 
down  at  15  Hz  to  reduce  high  frequency  noise.  The  resolution  is  1/4096  of 
these  values  or  0.09°,  0.012  ms-1  and  0.009  ns-1,  respectively. 

These  propellers  do  not  follow  a cosine  law  exactly,  but  compensation  for 
the  non-cosine  behaviour  can  and  must  be  made.  A detailed  discussion  is 
given  in  the  Appendix  along  with  a detailed  error  analysis.  Correction 
for  the  non-cosine  response  is  necessary  at  both  first  (important)  and 
second  order  (probably  not  essential  -see  Appendix).  The  propeller  with 
its  axis  nearly  horizontal  seems  to  follow  cosine  response  from  existing 
measurements  for  1 6 ] < 5-10°  and  measures 

Vj  = Q cos  6 + w sin  6 

where  5 is  the  small  tilt  angle  when  the  instrument  is  not  level,  with  6 
taken  positive  for  a tilt  down  into  the  wind  (see  Fig.  A.l)  and  Q is  the 
horizontal  wind  component  (U2  + v2)2.  The  propeller  with  its  axis  at  60° 
nominal  from  the  axis  of  the  vx  propeller  measures 

v2  = [Q  cos  (a  + 6)  + w sin  (a  + 6)  ] [1  - 0.328(6  - tan_1(w/Q) ) ] 

where  a is  the  angle  between  the  propeller  axes.  The  correction  factor 
1 - 0.328(6  - tan_1(w/Q))  comes  from  the  non-cosine  response.  There  is  some 
uncertainty  in  the  0.328  coefficient  (the  value  for  a = 60°)  and  its  varia- 
tion with  measured  values  of  a should  be  and  is  taken  into  account.  The 
deviation  from  cosine  response  for  the  measured  a (here  taken  to  be  60°) 
and  6 = 0 and  w = 0 is  included  in  the  calibration  for  v2  (see  Appendix): 

The  method  of  analysis  is  to  measure  and  72  (—  indicates  averages  over 

about  1 hour)  and  to  cake  v2-(72/71)v1  = -c  = aw2  + bw  and  solve  the 
quadratic  for  w.  The  effect  of  the  aw2  term  is  small  and  nearly  negli- 
gible as  shown  in  the  Appendix  but  is  retained.  The  various  errors  in 
calculating  w,  Q (and  hence  U = (Q2-v2)2)  are  discussed  in  the  Appendix. 
It  is  important  to  remember  that  any  errors  in  w and  6 calculated  from 
72/7i  will  propagate  into  the  calculated  q and  hence  u.  A summary  of 
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the  possible  errors  based  on  the  analysis  in  the  Appendix  is  given  in  Table  1. 
Careful  measurement  of  a is  clearly  essential.  Good  D.C.  stability  of  the 
amplifiers  and  offsets  (applied  to  make  the  outputs  - -5v  for  no  wind)  are 
also  important  and  we  have  been  very  careful  about  these  points  using  the 
best  low  temperature  coefficient  zener  diodes  for  reference  voltages.  The 
presence  of  some  Q in  the  tilted  propeller  also  helps  to  provide  a check 
of  proper  operation  since  low  pass  filtered  vx  and  v2  signals  should  and  do 
look  essentially  the  same.  There  is  no  double  check  on  the  wind  direction 
but  it  is  not  essential  to  measure  the  Reynolds  stress  (see  Appendix)  except 
to  indicate  when  we  are  in  the  wake  of  some  part  of  the  platform,  but  such 
effects  are  also  evident  from  the  vx  and  v2  measurements.  In  any  case  we 
have  only  had  to  replace  one  direction  potentiometer  when  a check  in  the 
laboratory  showed  it  was  not  working  well  over  its  full  range. 

d)  System  interconnections 

This  problem  is  often  not  discussed  but  in  fact  is  not  trivial.  The 
use  of  so-called  weather-proof  connectors  can  lead  to  grief  - often  con- 
siderable. In  our  system  we  use  underwater  connectors  (Mecca  connectors 
made  by  Teledyne  Corp.)  which  are  of  reasonable  cost;  other  similar 
connectors  are  also  available.  In  our  experience,  such  connectors  are 
essential  since  the  recording  package  part  of  the  system  goes  underwater 
at  times  (during  our  work  on  Sable  Island  and  on  the  Bedford  Institute 
stable  tower).  Thus,  we  recommend  the  use  of  underwater  connectors  in  all 
parts  of  the  system  exposed  to  the  marine  environment. 

e)  Sensor  response 

To  get  the  spectral  level  measurements  to  make  dissipation  estimates 
requires  a measurement  in  the  -5/3  range.  Because  of  the  response  limits 
of  the  sensors  we  must  work  at  the  low-frequency  end  of  the  range  to  avoid 
having  to  make  large  corrections  for  the  response  limits.  We  chose  to  use 
three  frequency  bands  centered  at  0.4,  0.8  and  1.6  Hz  as  a compromise. 

Having  three  bands  allows  a check  on  the  -5/3  slope  with  the  observations; 
it  also  allows  a band  to  be  dropped  if  its  n value  becomes  too  small 
(e.g.  the  0.4  Hz  band  when  Z/TJ  < h s)  or  when  the  response  correction 
becomes  rather  large  (e.g.  the  1.6  Hz  band  for  a D of  1 m and  winds 
less  than  10  ms-1)- 
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Table  1.  Summary  of 

possible  errors 

Cause 

Effect  on  |uw| 

Effect  on  CD 

Comments 

a error  of  t 1° 

± 5 % 

± 5 % 

We  believe  a error  is 

actually  within  ± 0 

and  effect  ± 2 - 3 %. 

Offset  error 

t 2 % at  5 ms-1 

± 1 % at  10  ms-1 

± 3 % at  5 ms-1 
± 2 % at  10  ms"1 

Actually,  offset  errors 
partly  cancel  rather 

than  add  and  are  about 

S these  values.* 

Calibration 

error  of  ± 2 % 

t 4 % 

± 1% 

Errors  in  B'(o) 

± 3 % 

± 3 % 

Fluctuations 

in  <5 

Negligible 

Negligible 

Non-cosine 
response  of  Vj 

t 1 % for  S = 0 

± 10  % for  5 = ± 10 

± 1 % for  S = 0 

± 10  % for  6 * ± 10 

For  \6\  < 2 - 3° 
similar  to  6 = 0. 

propeller  For  | 6 | =5  probably 

within  ± 5 %. 

If  everything  goes  wrong  including  all  the  data  having  5 equal  either 
+ or  - 5° , |uw|  may  be  in  error  by  about  ± 18*  and  Cp  by  ± 16%  but  with 
* varying  and  some  errors  perhaps  tending  to  cancel  rather  than  add,  we 
may  nope  to  get  average  Cp's  within  ± 10%  which  is  fairly  comparable  to 
any  other  system.  For  dissipation  estimates  only  calibration  and  offset 
instrument  errors  are  important.  For  a +2%  calibration  error  u;  is  in 
error  by  2.7%  and  Cp  by  -1.3%.  At  5 ms-1  the  offset  error  leads  to  a 0.3% 
error  in  u*  and  a 1.3%  error  in  CR.  Most  of  the  uncertainty  in  this  method 
comes  from  uncertainties  in  the  values  of  <,  K'  and  Z and  the  assumptions 
made  to  use  it. 

♦Note  that  the  quoted  values  require  d.c.  stability  of  1/2000  of  full  scale 
(50  ms"1  in  our  case)  for  2-3%  error  at  1/10  of  full  scale  (5  ms-1  in  our 
case) . 
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If  we  assume  that  the  sensors  behave  as  simple  RC  filters,  then  the  spectral 
values  from  two  different  bands  in  the  -5/3  range  can  be  used  to  determine 
their  response  characteristics  in  turbulence.  We  can  find  the  value  of  3 db 
down  frequency,  fo,  that  gives  the  same  value  of  c or  N>  from  the  two  bands, 
over  a range  of  wind  speeds.  If  the  response  can  be  represented  by  a distance 
constant  (0  = 2nfoD),  then  a plot  of  0 vs.  the  required  2n  fo  should  be  a 
straight  line  through  the  origin.  Figure  2 shows  such  a plot  based  on  hourly 
averages  of  the  0.8  Hz  and  1.6  Hz  filter  data  from  one  of  the  Gill  vi  pro- 
pellers on  the  Bedford  Institute  stable  platform.  Much  of  the  scatter  is 
probably  due  to  statistical  variations  from  the  -5/3  slope.  This  scatter 
increases  at  higher  values  of  0 because  the  corrections  become  relatively 
insensitive  to  fo  when  fo  is  large  compared  to  the  center  frequencies,  and 
therefore  it  takes  a larger  change  in  fo  to  compensate  for  the  same  relative 

difference  in  the  calculated  e's.  The  corrections  are  also  non-linear 

which  skews  the  scatter  to  high  values  of  2n  fo.  The  plots  suggest  that 

the  apparent  distance  constant  decreases  with  wind  speed.  Data  from  all 

three  filters  have  been  used  to  give  the  solid  line,  D = .56  0/(0-  1).  A 

fit  through  the  origin,  to  12  ms-1,  gives  a distance  constant  for  this 

anemometer  of  0.63  m.  Finding  a distance  constant  for  the  v_,  propeller 

is  more  complicated  because  we  essentially  measure  a combination  of  U and 

w.  One  must  use  observed  ratios  of  <t>  /$  and  4>  /t  to  estimate  the 

w u u w u 

departure  of  the  spectral  slope  from  -5/3.  The  value  calculated  for  the  V2 
propeller  is  .76  m of  air  flow  in  the  horizontal.  This  is  about  20%  larger 
than  the  vj  value,  perhaps  a reflection  of  the  deviation  from  cosine  response 
which  is  also  about  20*  for  a = 60°. 

We  have  found  similar  values  for  other  instruments,  however  one  of  these 
also  gave  values  of  1.0  m for  V]  and  1.2  m for  v?  for  some  early  data.  This 
was  a 'lit  puzzling,  but  it  turned  out  that  the  earlier  propellers  were 
nearly  twice  as  heavy  as  those  purchased  more  recently.  All  the  propellers 
we  have  used  are  the  7l,"  diameter  two-bladed  type  which  have  the  greatest 
strength.  The  quoted  distance  constant  is  0.8  ms*1  so  the  observed  performance 
of  the  lighter  propellers  in  turbulence  is  somewhat  better  than  claimed. 

As  noted  by  Hicks  (1972)  the  initial  response  is  a little  better  than 
expected  from  a single-pole  filter  behaviour  and  perhaps  this  behaviour 
explains  why  the  observed  distance  constant  is  better  than  that  predicted 
from  the  length  required  to  reach  (1-e*1)  of  full  speed  from  a zero  start. 
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It  is  no  longer  possible  to  obtain  the  type  of  propeller  we  have  been  using 
- a new  four-bladed  version  with  slightly  different  calibration  is  now  being 
made.  We  have  not  established  its  response  as  yet,  but  judging  from  the 
weight  it  should  be  similar  to  the  better  old  type  of  propeller. 

Clearly,  if  response  corrections  are  needed,  a way  to  establish  the  response 
under  the  conditions  of  operation  is  necessary.  With  the  response  corrections 
the  drag  coefficients  calculated  from  the  stresses  estimated  from  the  three 
bands  generally  agree  with  the  mean  of  the  three  values  within  about  5*  for 
winds  from  6-20  ms-1  and  the  ratios  between  stresses  from  pairs  of  bands 
are  independent  of  wind  speed  on  average. 

A similar  approach  was  used  to  determine  the  response  of  the  microbead 
thermistors  in  their  protective  enclosure.  We  find  that  they  effectively 
have  a distance  constant  of  1 metre.  Thus,  their  response  is  being  limited 
by  the  enclosure  which  we  feel  is  necessary  for  protection  from  rain  and 
spray  - the  radiation  shielding  could  probably  be  done  less  restrictively . 
However,  the  observed  response  is  adequate  for  the  measurements  we  want  to 
make,  although  we  had  hoped  it  would  be  somewhat  better. 

IV.  RECORDING  SYSTEMS 

The  data  required  for  Reynolds  flux  and  dissipation  estimates  of  fluxes  are 
recorded  by  two  separate  systems. 

a)  Reynolds  flux  system 

The  Reynolds  flux  cospectra  span  three  to  four  decades  and  require 
recording  a great  deal  of  data.  The  amount  of  data  recording  can  be  greatly 
reduced  if  we  do  not  sample  the  higher  frequencies  all  the  time.  Sub- 
sampling will  increase  the  statistical  uncertainty  a little  but  makes  it 
possible  to  obtain  many  more  estimates  and  therefore  improves  the  overall 
statistics.  Note  that  statistical  uncertainty  in  the  low  frequencies  is  an 
important  source  of  the  variaoility  of  the  bulk  aerodynamic  coefficients. 

To  establish  their  behaviour  with  various  parameters  (e.g.  wind  speed  or 
stability)  requires  a large  number  of  estimates.  For  remote  operation,  sub- 
sampling is  required  because  of  the  data  storage  limitations  of  low  power 
recording  systems. 
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Fig.  3.  Block  diagram  of  data  flow  in  the  Reynolds  flux 
recording  package. 


A block  diagram  of  the  data  flow  in  this  recording  system  is  shown  in  Fig.  3. 
Up  to  six  channels  may  be  sampled.  Three  of  these  are  from  the  twin 
propeller  vane  anemometer;  the  other  three  are  scalar  signals  - the  micro- 
bead thermistor  and  two  of  the  various  humidity  sensors  we  have  tried. 

During  our  present  work  on  the  station  P weathership  CGSS  Quadra,  we  have 
used  two  microbeads.  In  our  planned  work  during  JASIN  we  shall  record  a 
microbead,  a-Lyman  humidity  and  dew  point  temperature.  By  recording 
these  two  humidity  sianals  we  car  calibrate  the  a-Lyman  against  the  dew 
point  at  higher  frequencies  than  we  can  with  the  dissipation  system. 

In  the  Reynolds  flux  system  each  signal  from  a transducer  preamplifier 
goes  to  a pair  of  circuits:  one  is  a low-pass  filter  with  a one-second 
time  constant  (3  db  down  at  0.16  Hz);  the  other  is  a prewhitening  circuit 
to  emphasize  the  higher  frequencies.  The  low-pass  filter  is  an  RC  filter 
with  a micropower  operational  amplifier  follower  on  its  output  (gain 
1.000).  This  signal  is  used  for  a slow  sampling  rate  sample.  The  pre- 
whitening circuit  acts  like  a time  differentiator  at  low  frequencies, 
reaches  maximum  response  at  1.5  Hz  and  then  rolls  off  like  an  RC  low-pass 
filter  at  higher  frequencies.  It  also  uses  a micropower  operational 
amplifier  as  the  active  element.  This  signal  is  used  for  a fast  sampling 
rate  sample. 
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If  no  low-pass  filtering  were  done  and  the  signals  were  sampled  at  a fairly 
slow  rate,  the  higher  frequencies  would  alias  back  and  variance  and  covariance 
should  be  preserved.  However,  the  spectra  and  cospectra  could  be  very  dis- 
torted. If  unusual  results  were  obtained  (e.g.  a dramatic  jump  in  Cq  at 
high  wind  speed)  it  would  be  difficult  to  tell  what  the  cause  might  be. 

Also,  if  the  higher  frequencies  require  instrument  response  corrections 
they  cannot  be  applied  unless  the  spectral  shape  is  assumed.  We  chose 
to  low-pass  filter.  With  our  usual  slow  sampling  rate  (1/3  Hz,  i.e.,  one 
sample  every  three  seconds)  the  little  bit  of  aliasing  occurring  after 
filtering  should  nearly  compensate  for  the  filter  loss.  The  approach  of 
allowing  full  aliasing  and  using  the  high  frequency  subsampling  to  tell 
what  aliasing  is  occurring  and  hence  what  instrument  corrections  are  needed 
should  work  but  is  somewhat  more  complicated. 

By  using  the  prewhitening  circuit  we  eliminate  any  spectral  distortion  due 
to  frequencies  below  those  sampled  (which,  because  of  the  greater  amplitudes 
at  low  frequencies,  have  an  effect  similar  to  a large  trend).  The  signal 
level  is  also  increased,  which  keeps  the  digitization  noise  at  a very  low 
level.  The  gain  is  chosen  so  that  the  mean  square  output  is  about  1/50  of 
full  scale  for  the  largest  expected  signals  allowing  fairly  large  peak 
factors  and  considerable  error  in  our  estimates  of  the  required  gain  (for  a 
sine  wave  the  peak-to-peak  value  would  be  1/5  of  full  scale).  In  practice 
overloading  has  not  occurred  so  far  in  our  observations.  Our  fast  sampling 
rate  is  3 Hz  (limited  by  the  digital  recorders).  By  allowing  the  response 
to  be  flat  at  1.5  Hz  (the  Nyquist  frequency)  we  do  allow  the  high  frequencies 
which  we  cannot  sampie  to  alias  back  so  that  we  do  not  lose  their  contribu- 
tion. Of  rnur<;p  the  frequencies  above  1.5  Hz  are  somewhat  undercorrected 
for  inst/ument  response.  Except  for  4>w , the  aliasing  effects  are  small. 
Undercorrection  for  instrument  response  should  be  negligible  except  perhaps 
for  o . Since  our  primary  concern  is  the  fluxes  these  limitations  are  not 
serious. 

Sample  and  hold  circuits  are  used  to  avoid  the  need  to  correct  for  time 
skewness  in  the  sampling.  The  multiplexer  allows  the  A/D  (analog-to-digital 
converter)  to  sample  the  fast  or  slow  channels  in  sequence  (up  to  a maximum 
of  6). 
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The  A/D  converter  is  a 12-bit  low  power  unit  (Analog  Devices).  With  12-bit 
accuracy  the  signal  to  noise  ratio  is  potentially  about  83  db.  In  our 
experience,  digitization  noise  is  the  limiting  factor;  with  no  signal  (no 
propellers  and  dummy  resistors  for  temperature  and  humidity  sensors)  the 
change  in  the  A/D  value  is  at  most  ± 1 bit.  Now  the  rms  fluctuations  for 
wind  are  about  1/10  of  the  mean,  or  about  20  db  down  at  full  scale  (50  ms-1), 
and  40  db  down  at  1/10  full  scale  (5  ms-1).  For  the  slow  samples  the  vari- 
tion  over  the  frequency  range  sampled  is  perhaps  another  20  db  so  even  the 
lowest  spectral  values  have  about  20  db  signal-to-noise  ratio.  For  the 
fast  samples  the  signal  is  about  17  db  below  maximum  at  50  ms-1  and  44  db 
at  5 ms-1.  The  spectra  are  nearly  flat  so  there  should  be  no  problem. 

For  temperature,  the  full  scale  is  about  20  C°;  with  an  air-sea  temperature 
difference  of  1°  the  rms  fluctuations  are  about  0.1°,  or  about  46  db  down 

for  the  slow  samples  and  38  db  down  for  the  fast  samples  if  l)  = 5 ms-1. 

In  the  scheme  we  have  usually  employed,  the  prewhitened  signals  are  sampled 
at  3 Hz  for  128  samples  (42.7  sec)  and  then  the  low-pass  filtered  signals 
are  sampled  at  1/3  Hz  for  256  samples  (12.8  minutes).  This  sequence  is 
repeated  four  times  to  form  one  data  group  to  give  one  Reynolds  flux  estimate 
(over  a period  of  54  minutes).  A time  word  (the  number  of  0.1-hour  intervals 
since  start-up)  is  written  at  the  end  of  each  group.  As  it  is  also  12  bits 
it  folds  over  after  17  days,  but  the  folds  are  easy  to  allow  for.  A check 
on  the  contributions  from  lower  frequencies  can  also  be  made  from  the  means 
over  the  12.8  minute  slow  samples.  The  fast  sample  rate  is  fixed  but  the 

slow  sample  rate  may  oe  varied  from  1/3  sec  to  33  sec  in  1/3-sec  increments. 

The  number  of  samples  of  each  type  may  be  varied  from  2‘4  (16)  to  213  (8192) 
and  the  number  of  fast  and  slow  sample  pairs  in  a group  from  1 to  9.  The 
time  interval  between  collecting  groups  may  be  varied  from  0.1  to  99.9  hours 
in  0.1-hour  increments.  Usually  we  have  used  3.0-hour  intervals. 

The  digital  cassette  recorders  we  have  used  are  those  made  by  Memodyne. 

Their  specifications  allow  them  to  be  used  at  up  to  100  stpps  (and  bits) 
per  second.  To  get  our  3 Hz  sample  rate  wo  run  a pair  of  recorders  at  160 
steps/sec  to  get  three  channels  per  recorder.  At  this  stepping  rate  they 
are  somewhat  sensitive  and  one  must  pretest  the  cassettes  to  he  sure  they 
will  work  well  in  a particular  machine.  In  retrospect,  it  might  have  been 
better  either  to  have  accepted  a 2 Hz  sampling  rate  and  a little  more  aliasing 
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distortion  or  to  have  run  three  recorders  with  two  channels  each.  We  record 
the  12-bit  words  as  a 2-bit  gap  (no  flux  changes)  and  the  12-bit  word.  A 
mark  character  (2-bit  gap  + 1 bit)  separates  each  scan  (sample)  of  three 
channels.  These  choices  were  made  to  avoid  losing  track  of  the  words  even 
if  a bit  were  lost.  In  practice,  we  have  had  more  trouble  with  gap  recog- 
nition than  losing  bits  and  could  probably  have  omitted  the  gaps  completely, 
but  we  didn't  know  what  the  limits  were  when  we  started. 

Originally  we  started  with  standard  digital  cassettes  (equivalent  to  a C60 
audio  cassette  in  length).  Later  we  switched  to  C90  type  cassettes  (Informa- 
tion Terminals)  to  increase  the  record  length.  C120  cassettes  have  also 
been  used  by  others  but  with  our  high  stepping  rate  might  not  be  satis- 
factory - we  have  not  made  a systematic  check  on  such  cassettes. 

With  a C60  cassette  we  have  room  for  30  groups  of  four  each  of  128  fast  and 
256  slow  scans;  with  the  C90  cassettes  we  get  44  such  groups  on  one  tape. 

With  a 3-hour  interval,  this  allows  5.5  days  of  data  on  a C90  cassette.  To 
extend  the  data  length  further  we  use  three  pairs  of  recorders  which  allows 
16.5  days  of  recording. 

Even  this  length  is  somewhat  short.  More  recorders,  a different  type  of 
recorder  or  a longer  interval  between  recordings  provide  alternatives. 

Sampling  at  greater  than  3-hour  intervals  is  not  attractive  since  one  would 
get  only  a very  small  number  of  samples  during  a particular  high-wind  speed 
event.  (Sampling  more  often  is  probably  unnecessary  as  more  frequent  samples 
are  not  completely  statistically  independent.)  Recognition  of  the  fact  that 
many  of  the  samples  will  not  be  of  interest  allows  one  to  extend  the  total 
sampling  length.  Low  winds  ( ^5  ms-1)  are  of  little  interest.  Even  winds 
< 10  ms*1  soon  become  of  little  interest.  A circuit  was  added  to  average 
the  wind  speed  over  0.1-hour  (six  minute)  intervals,  the  basic  interval 
rate  of  the  recording  system.  When  a record  was  to  be  taken,  it  would  not 
be  allowed  if  the  mean  speed  was  below  a set  limit  in  the  preceding  six  minutes. 
The  settings  were  0-9,  allowing  recordings  for  all  winds  at  setting  0 and 
6 to  22  ms*1  in  2 ms*1 -increments  for  settings  1 to  9.  Thus,  with  a setting 
of  3 (10  ms*1  lower  limit)  as  usually  used  on  the  Bedford  stable  platform 
experiment,  records  would  be  taken  only  about  1/3  of  the  time  on  average, 
giving  recordings  for  up  to  about  45  days.  With  an  intended  service  interval 


of  30  days,  this  scheme  seemed  to  provide  sufficient  length  to  allow  for 


bad  weather  delaying  servicing  and  proved  adequate  in  practice. 

The  prewhitened  signals  must  be  corrected  for  instrument  response.  In  the 
case  of  the  v:  and  v2  signals  these  corrections  must  be  applied  before  u 
and  w can  be  calculated  because  of  the  need  for  a non-linear  correction 
(see  Appendix).  The  Vj  and  v2  signals  are  Fourier  transformed,  corrected 
for  the  circuit  and  instrument  response  and  inverse  transformed  using  the 
first  slow  sample  value  following  the  fast  sample  to  give  the  offset  needed 
for  Vj  in  calculating  w.  u and  w are  then  calculated  and  Fourier  trans- 
formed. To  calculate  uw,  other  covariances  involving  w and  spectra,  the 
slow  sample  values  are  used  for  frequencies  up  to  their  Nyquist  frequency 
(1/6  Hz)  and  the  fast  sample  spectral  and  cospectral  values  at  frequencies 
above  1/6  Hz  are  used  for  the  high  frequencies.  Originally  we  planned  to 
'correct1  the  fast  sample  values  by  comparing  spectral  and  cospectral  values 
in  the  overlap  region  1/43  Hz  to  1/6  Hz,  but  because  of  statistical  variation 
in  the  lower  frequencies  of  the  fast  sample  results  this  procedure  seems  to 
be  less  satisfactory  than  just  using  the  observed  fast  sample  values  for 
frequencies  above  1/6  Hz.  More  sampling  to  reduce  the  variations  above  1/6  Hz 
should  be  done  if  possible,  although  we  couldn't  during  the  Bedford  stable 
platform  experiments. 

b)  Dissipation  system 

This  package  provides  continuous  recording  of  the  variables,  their  low- 
frequency  variations,  and  spectral  levels  from  which  dissipation  estimates 
may  be  made.  The  data  flow  for  a channel  for  which  both  low-frequency  and 
spectral  band  information  is  desired  is  shown  in  Figure  4. 
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Fig.  4 Data  flow  in  the  dissipation  recording  system. 

Spectral  band  data  are  obtained  for  both  the  Vj  and  v?  propellers.  The  v2 
propeller  data  provide  a check  and  are  collected  to  allow  dissipation  esti- 
mates to  be  made  if  the  Vj  channel  fails  before  the  v2  channel.  Such  esti- 
mates being  a mixture  of  4^,  4>w  and  <j>uw  would  have  more  uncertainty;  how- 
ever, so  far,  when  we  have  data,  both  signals  are  present  and  only  the 
band  pass  values  are  used.  Spectral  band  data  are  also  collected  from  the 
microbead  and  have  been  attempted  (unfortunately  with  little  success)  from 
the  various  humidity  sensors.  Hopefully,  we  shall  get  some  useful  data 
with  the  a-Lyman  system  during  JASIN. 

The  low  pass  filters  are  25-second  single-pole  RC  circuits  followed  by  a 
micropower  operational  amplifier  in  a follower  mode  (gain  1.000).  Channels 
for  which  spectral  data  are  not  required  have  the  low-pass  filter  treatment 
only,  e.g.  sea  temperature,  wind  direction,  filtered  Brady  array  (when  used), 
air  temperature  from  the  rod  thermistor,  air  temperature  from  the  dew  point 
system  (when  used)  and  dew  point  temperature  (when  used).  A maximum  of  20 
low-pass  filter  channels  can  be  recorded  but  up  to  now  we  have  only  recorded 


28 


nine  data  channels  plus  a zero  reference  to  check  for  drift.  In  our  earlier 
work  the  zero  reference  was  a grounded  low-pass  filter  circuit  (no  drift 
within  the  A/D  resolution  was  experienced).  When  we  added  channels  and  used 
all  of  the  nine  low-p  Filter  circuits  initially  installed  the  zero  refer- 
ence became  a grounded  multiplexer  channel. 

The  prewhitener  circuit  acts  as  a differentiator  at  low  frequencies,  reaches 
peak  gain  at  about  10-  15  Hz  and  rolls  off  as  a single-pole  RC  filter  at 
higher  frequencies.  A micropower  operational  amplifier  is  the  active  element. 
It  makes  the  input  spectrum  to  the  band  pass  filters  nearly  white.  For  a 
-5/3  input  spectrum  the  output  spectrum  is  +1/3  in  the  band  pass  filter 
range.  The  gain  was  chosen  to  make  the  mean  square  output  about  1/50  of 
full  scale  for  maximum  input  (e.g.  50  ms-1  for  wind,  aT  = 5 C°  with 
U = 25  ms"1  for  temperature).  The  band  pass  filters  are  an  RC  double-pole 
high-pass  stage  followed  by  a double-pole  low-pass  stage  using  a dual  micro- 
power operational  amplifier  as  the  active  element.  Care  was  taken  to  make 
the  final  DC  output  very  small  and  stable.  Again,  the  gains  were  chosen 
assuming  a -5/3  original  input  spectrum  to  make  the  mean  square  output 
1/50  full  scale  at  maximum  input  to  allow  fairly  large  peak  factors  and 
leave  leeway  for  errors  in  the  design  criteria  (for  a sine  wave  the  peak- 
to-peak  value  would  be  1/5  of  full  scale). 

Over  most  of  the  interval  between  recordings  each  band  pass  channel  is 

sampled  at  20  Hz,  the  output  digitally  squared  and  summed.  (In  a second 

package  using  a microprocessor,  the  sample  rate  is  10  Hz  because  squaring 

by  the  microprocessor  is  slower  than  the  hard-wired  version  of  the  original 

package).  Analog  squaring  of  the  signals  was  considered  briefly,  but  was 

quickly  discarded  as  having  inadequate  dynamic  range.  The  accurate  range 

of  analog  squaring  is  about  20  db  or  so.  However,  for  a wind  velocity 

signal  with  Cn,  a ,,/u*,  and  hence  o / 0 independent  of  U,  the  output  power 
u w u —8/3 

of  a fixed  frequency  band  pass  filter  is  proportional  to  U . For  a 

change  in  U from  5 to  50  ms'1  the  power  output  changes  by  about  500,  or  27  db. 

Since  we  must  allow  about  20  db  from  full  scale  for  peaking  and  possible 

underestimates  of  Cn  at  high  U,  analog  squaring  is  much  too  inaccurate. 

2/3 

For  scalars  the  power  output  is  proportional  to  U (a>) 
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In  calculating  the  squares,  one  has  effectively  an  11-bit  number  since  the 
sign  does  not  matter.  We  decided  to  retain  the  high  order  20  bits  of  the 
square,  the  lower  2 bits  probably  being  mainly  noise.  These  20  bits  were 
accumulated  in  a 32-bit  register.  With  4096  squares  (^  205  seconds),  full- 
scale  readings  would  fill  the  32-bit  register  but  since  the  mean  square 
value  should  be  about  1/50  of  full  scale  the  register  will  have  zeros  in 
the  first  5 or  6 bits  in  fact. 

To  record  the  values  in  the  registers  a floating  point  format  was  used. 

Shifting  along  the  register  for  up  to  15  shifts  was  used  until  a one  bit 
was  found.  The  number  of  shifts  was  recorded  in  the  first  4 bits  and  the 
8 bits  in  the  register  after  shifting  was  recorded  as  the  last  8 bits. 

This  scheme  does  not  use  the  full  dynamic  range  of  the  system,  in  fact,  and 
also  does  not  give  a full  12  bits  of  information  except  for  small  values. 

If  there  are  less  than  15  shifts  the  first  bit  in  the  remaining  eight  is 
always  a one.  However,  the  scheme  used  does  allow  adequate  dynamic  range 
and  ensures  that  non-zero  bits  are  signal,  not  noise.  A shift  of  up  to 
25  places  (requiring  a 5-bit  exponent)  could  be  used  if  more  dynamic  range 
is  needed,  but  the  resolution  of  the  mantissa  would  only  be  1/64  unless 
the  first  bit  was  dropped  except  when  25  shifts  occurred.  The  dynamic 

range  allowed  is  224  or  72  db,  less  than  the  noise  level  of  83  db  which  will 

be  further  reduced  by  summing  many  values. 

Suppose  F = 2.5  ms-1  or  1/20  of  full  scale;  then  the  power  output  of  a 
band  pass  filter  of  a velocity  signal  is  about  (20)  °^,  or  35  db  down  from  the 
value  expected  at  50  ms-1  which  is  17  db  below  full  scale,  for  a sum  of 

4096  squares;  the  52  db  down  total  is  still  20  db  above  the  zero  value.  If 

a larger  or  smaller  number  of  squares  are  summed,  there  is  still  plenty  of 
range.  For  example,  for  1024  iterations  the  signals  at  U = 2.5  ms'1  are  still 
14  db  above  the  zero  value,  i.e.  the  resolution  is  still  1/25  of  the  reading. 
For  temperature  with  U = 1.25  ms-1  and  aT  = 0.25  C°  (both  1/20  of  design) 
one  is  also  20  db  above  the  minimum  value,  for  a sum  of  4096  squares. 

Thus,  the  scheme  we  have  used  is  adequate  since  the  minimum  values  used  as 
examples  are  below  those  used  in  practice,  and  when  one  gets  near  zero 
values  one  is  sure  the  fluxes  are  small  (as  confirmed  by  the  small  directly 
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recorded  values  of  U and  .\T).  For  very  light  winds  ( «•  1-2  ms'1)  brush 
noise  from  the  tacogenerators  may  give  non-zero  band  pass  outputs  for  Vi 
and  v2»  but  these  are  obvious  and  brush  noise  is  at  too  high  a frequency 
to  cause  measurable  effects  at  wind  speeds  of  interest. 

The  low-pass  filter  outputs  and  the  sums  of  squares  of  band  pass  filter 
outputs  (up  to  20  channels)  are  recorded  at  intervals.  One  to  20  of 
each  type  are  selectable  but  9 low-pass  and  15  band  pass  filter  channels 
plus  a zero  reference  for  each  are  the  maximum  used  so  far.  After  each 
scan,  a 4-bit  identification  word  and  the  lower  order  8 bits  of  the  scan 
count  are  recorded.  The  high  order  12  bits  of  the  scan  count  is  also 
recorded  after  each  block  of  scans  (maximum  number  before  folding  is  2:’°). 

The  number  of  scans  per  block  may  be  selected  from  1 to  99  and  the  number 
of  blocks  per  tape  from  1 to  999. 

The  sample  interval  may  be  varied  from  about  6 sec  (minimum  to  allow 
recording  and  bookkeeping)  to  999  sec.  The  number  of  iterations  (number  of 
squares  in  the  sum)  is  chosen  to  use  most  of  the  time  between  recordings. 
Maximum  allowed  is  9999  at  20  Hz  (8.3  minutes)  in  the  first  package  and 
4090  at  10  Hz  (6.8  minutes)  in  the  second  package  built.  Our  basic  design 
time  interval  was  4 minutes,  allowing  up  to  4600  squares  to  be  summed  (or 
2300  in  the  microprocessor  controlled  package) . We  have  used  2-  or  4-minute 
intervals  in  earlier  experiments.  Presently  on  the  weathership  program,  we 
are  using  5-minute  intervals  allowing  28  days  of  data  on  one  tape.  We  use 
two  recorders  in  sequence,  allowing  56  days  total  which  is  a bit  longer 
than  one  patrol  of  48  days  including  time  to  and  from  station. 

c)  Circuit  response  measurements 

The  detailed  response  of  many  circuits  must  be  measured  in  the  system  - 
in  particular,  the  6 prewhitener  circuits  in  the  Reynolds  flux  package  and 
the  15  band  pass  filters,  including  the  effect  of  the  preceding  prewhi teners , 
in  the  dissipation  package.  Measuring  the  response  by  analog  means  is  both 
tedious  and  awkward  because  of  the  frequencies  involved.  A good  chart 
recorder  is  about  the  only  means  and  amplitude  ratio  accuracy  of  about  2'1 
is  about  the  best  one  can  do.  This  method  was  used  at  a few  points  to  check 
the  digital  method  actually  used. 
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The  method  used  was  to  generate  a white  spectrum  (sum  of  sine  waves  of 
equal  amplitude)  over  the  frequency  range  of  interest  on  our  PDP-12  computer 
and  generate  an  analog  signal  using  a D/A  converter.  The  D/A  holds  the  out- 
put until  updated  and  thus  has  a staircase  character,  but  by  going  to  some- 
what higher  frequency  than  finally  needed  and  low-pass  filtering,  a white 
spectrum  with  a smooth  roll-off  is  produced.  This  spectrum  is  passed  through 
the  circuits  and  both  the  input  and  output  are  digitized  on  the  PDP-12.  Fast 
Fourier  transforms  are  done  on  the  UBC  Computing  Centre  IBM  370  from  which  the 
transfer  function  including  phase  nay  be  computed  accurately.  As  the  PDP-12 
has  a 10-bit  A/D  the  dynamic  range  is  60  - 70  db  which  is  adequate. 

Originally  we  tried  generating  white  noise  using  random  numbers,  but  the 
white  spectrum  with  the  generation  frequencies  matched  to  the  sampled 
frequencies  gives  more  coherence  between  input  and  output  when  the  transfer 
function  is  small. 

For  the  prewhitening  circuits  of  the  Reynolds  flux  package  the  transfer 
function  is  quite  small  at  low  frequencies  but  needs  to  be  known  accurately. 

A straight  line  was  fitted  in  this  region  to  smooth  out  the  small  random 
errors.  For  the  band  pass  filters  of  the  dissipation  package  such  measures 
are  not  needed.  However,  because  the  bandwidth  is  fairly  large  it  is 
necessary  to  allow  for  the  spectral  shape  of  the  actual  turbulence  data. 

In  calculating  the  relation  between  the  input  power  and  the  final  output,  a 
-5/3  spectrum  was  assumed. 

d)  Power  requirements. 

Digital  recording  offers  great  advantages  for  remote  operation  in  terms 
of  wide  dynamic  range  and  low  power.  (Although  the  low  power  digital  tape 
units  require  some  compromise,  e.g.  subsampling  the  higher  frequencies  in 
the  Reynolds  flux  package.)  The  overall  power  requirement  for  the  trans- 
ducers, the  flux  package  and  the  dissipation  package  is  about  50  mA  at  36  v 
nominal  (28  - 40  v range).  The  36  v was  chosen  to  be  compatible  with  the 
power  supply  on  the  Bedford  Institute  of  Oceanography  stable  platform;  the 
system  could  be  adapted  to  run  on  24  v fairly  easily  to  reduce  the  number  of 
batteries  required.  Negative  supply  voltages  needed  for  the  analog  circuits 
are  generated  internally.  Substitution  of  the  microprocessor  dissipation 
package  increases  the  current  drain  to  about  100  mA,  so  the  greater  flexibility 
of  using  a microprocessor  does  require  extra  power. 


V.  IN-FIELD  READER 

A serious  problem  with  all  digital  recording  instruments  is  how  to  check 
them  for  proper  operation  in  the  field.  Monitoring  connectors  on  the 
recording  packages  were  included  so  the  preamplified  transducer  signals 
could  be  checked.  One  can  also  see  if  the  tape  recorders  are  advancing 
correctly  but  these  checks  are  insufficient  and  it  is  not  satisfactory  to 
wait  until  the  cassette  can  be  translated  to  standard  9-track  tape  on 
our  PDP-12  system  in  the  laboratory.  Thus,  to  provide  a check  in  the 
field,  a reader  unit  was  constructed  to  decode  the  output  from  a Memodyne 
reader  unit  and  provide  signals  for  a six-channel  chart  recorder  (a  Brush 
260  in  our  case) . 

For  the  dissipation  package  records  the  unit  can  be  set  to  start  on  any 
channel  and  decode  six  sequential  channels.  A 10-bit  D/A  was  used  because 
a 12-bit  unit  was  not  readily  available  when  the  unit  was  built,  but  it  is 
adequate  for  chart  monitoring.  The  low-pass  channels  are  put  out  directly. 

For  the  band  pass  channels  the  first  4 bits  (the  exponent)  are  inverted  so 
that  increasing  power  gives  increasing  output,  the  first  bit  of  the  remaining 
8 is  dropped  (unless  the  exponent  indicates  15  shifts)  and  the  others  are 
shifted.  The  output,  then,  is  approximately  log2  (of  band  pass  output). 

Actually,  it  is  a set  of  sixteen  linear  segments  approximating  a logarithmic 
curve.  The  in-field  reader  proved  invaluable  in  checking  for  proper  opera- 
tion and  gives  an  excellent  indication  of  it.  For  example,  when  working 
properly  the  Vi  and  v2  low-pass  signals  should  and  do  look  very  similar, 
the  microbead  and  rod  thermistor  air  temperatures'  low-pass  signals  look 
similar.  Likewise,  the  band  pass  outputs  for  a given  transducer  all 
follow  one  another.  The  reader  will  also  put  out  the  three  data  channels 
plus  time  channel  from  a Reynolds  flux  package  tape. 

VI.  FIELD  MEASUREMENTS  AND  SOME  PRELIMINARY  RESULTS 

After  local  field  testing  our  first  field  expedition  with  the  equipment 

was  to  Sable  Island,  Nova  Scotia,  for  an  intercomparison  experiment  with 

the  air-sea  interaction  group  of  the  Bedford  Institute  of  Oceanography  in 

September  and  October  1975.  Our  preamplified  transducer  signals  were  also  * 

recorded  on  the  BIO  recording  system  and  initial  analysis  of  these  records 

was  performed  at  BIO  for  comparison  with  their  results.  ] 
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The  shortcomings  of  the  humidity  sensors  we  were  trying  which  this  experi- 
ment revealed,  have  already  been  mentioned  in  section  III  b.  A need  to  im- 
prove our  temperature  measurements  was  also  indicated,  as  discussed  in 
section  111  a. 

The  results  of  this  experiment  were  written  up  jointly  as  a BIO  technical 
report  (Smith,  et  al.  1976).  The  main  results  are  that  all  the  wind 
sensors  (Sonic  and  thrust  from  BIO  and  Gill  from  10UBC)  gave  comparable 
results  (our  GILL  anemometer  system  required  response  corrections,  as 
expected)  and  that  the  indicated  values  of  Cp  are  similar  to  those  reported 
by  Smith  and  Banke  (1975)  mainly  based  on  offshore  tower  data.  Detailed 
response  corrections  to  the  GILL  were  not  made  but  were  allowed  for  by 
comparison  with  the  BIO  results  and  seemed  reasonable  if  the  distance 
constant  were  about  1 m.  At  that  time  we  had  not  established  the  exact 
distance  constants  as  discussed  in  section  III  e.  However,  the  covariances 
at  a given  height  come  from  a fixed  range  of  length  scales  so  if  the  response 
limit  occurs  as  a distance  constant,  the  correction  is  independent  of  0.  Our 
own  analysis  of  these  results  is  still  in  progress  but  it  appears  that  with 
detailed  corrections  both  Reynolds  flux  and  dissipation  estimates  give 
similar  results  to  those  presented  in  Smith  et  al.  (1976). 

The  next  field  experiment  was  also  conducted  with  the  air-sea  group  at  BIO 
on  their  stable  platform  - about  10  miles  offshore  of  Halifax,  Nova  Scotia. 

Again,  analysis  on  their  system  gives  comparable  results  with  smaller 
response  corrections  to  the  GILL  system  indicated  (probably  because  of  the 
greater  height  and  also  because  of  the  better-response  newer  propellers). 

Our  analysis  is  giving  similar  results  with  Cp  values  generally  agreeing 
wi th  the  Smith  and  Banke  (1975)  curve  showing  an  apparently  fairly  slow 
linear  increase  of  Cp  with  0.  Measurements  were  attempted  from  September 
1976  through  April  1977  and  a great  deal  of  data  were  obtained.  We  are 
particularly  pleased  with  the  amount  of  data  obtained  above  15  ms-1.  Table 
2 shows  the  number  of  Reynolds  flux  estimates  and  number  of  hours  of  dissi- 
pation data  obtained  for  unlimited  fetch  conditions  (*2-  1 hour  of  dissipa- 
tion data  is  sufficient  for  a flux  estimate).  There  are  also  data  for 
limited  fetch  which  we  shall  examine.  We  also  plan  to  look  at  the  time  variation 
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of  Cq  with  storm  front  passages,  13  of  which  we  have  identified. 


Table  2.  Unlimited  fetch  data  from  the  Bedford  Institute  stable  platform. 


wind  speed  (ins'1) 

00 

1 

<sO 

9 - 11 

12  - 14 

15  - 17 

18  - 20 

# of  Re  flux  records 

67 

79 

82 

32 

13 

# of  hrs  of  diss.  data 

165 

175 

77 

53 
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Our  next  field  experiment  started  in  July  1977  with  our  system  being 
installed  on  board  CCGS  Quadra,  one  of  the  two  ships  which  occupy  ocean 
weather  station  PAPA.  The  first  patrol  was  during  MILE  (Mixed  layer  experi- 
ment) and  in  spite  of  unintended  shut-down  problems  caused  by  power  line 
transients,  a good  deal  of  data  were  obtained  thanks  to  the  efforts  of 
Dr  M.  Miyake  and  his  staff  who  were  on  board  running  another  program.  A 
good  deal  of  work  had  to  be  done  to  reduce  radio  transmission  interference 
which  was  increased  by  adding  the  a-Lyman  and  dew  point  systems  and  their 
cables.  After  the  second  patrol  the  dew  point  and  a-Lyman  systems  were 
removed  as  it  became  obvious  that  few  useful  data  were  being  obtained  by 
them.  On  the  third  patrol  which  was  recently  completed,  the  system  ran 
throughout  although  both  microbeads  were  broken  after  a few  days.  The 
last  patrol  on  which  we  shall  operate  starts  in  late  March  1978  and  ends 
in  early  May.  Quite  a few  observations  for  winds  up  to  20  ms'1  have  been 
obtained  with  a few  occasions  of  winds  to  25  ms-1.  An  exact  list  including 
corrections  for  ship's  velocity  has  not  yet  been  compiled. 

Our  next  field  operation  (and  probably  the  last,  at  least  for  some  time) 
will  be  on  the  Meteor  during  JASIN.  Here,  one  of  us  (WGL)  will  be  on 
board  and  another  attempt  can  be  made  to  get  moisture  flux  measurements 
using  the  a-Lyman  humidiometer  for  fluctuations  with  the  dew  point  system 
to  provide  in  situ  calibration. 
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Appendix.  Analysis  of  the  Gill  twin  propeller  vane  anemometer. 


Fig.  A.l.  Sketch  of  propeller  system  to  define  angles; 
6 is  tilt  of  Vj  propeller  (tilt  down  taken  to  be  positive), 
a is  angle  between  propeller  axes,  nominally  60°,  e is  angle 
of  attack. 


The  propellers  are  vane-mounted  which  keeps  their  axes  and  the  instan- 
taneous wind  vector  in  the  same  plane.  If  the  propellers  were  cosine 
devices  so  that  the  axial  velocity  component  v..  = total  velocity  times 

cosine  of  angle  of  attack,  then 
vt  = Q cos<S  + w sins 

v'  = Q COs(a+6 ) + w sin(a+^) 

where  Q = Q + q,  mean  plus  fluctuation  of  the  horizontal  velocity  component 
and  w is  the  vertical  velocity  component.  Now  Q = [(U  + u)2  + v2]’*  where 


U = U + u is  the  downstream  component,  v is  the  cross-stream  component 


and  u = v = 0.  Thus  expanding  Q = U + u+^+  higher  order  terms 


oy  = ( v 2 ) i = U/ 10  so  the  v2  term  is  quite  small  and  furthermore  will  be 


almost  uncorrelated  wi  th  u or  w and  may  usually  be  neglected.  In 
practice,  we  resolve  Q into  l)  + u and  v using  the  wind  direction 
indicated  by  the  vane,  but  if  the  vane  circuitry  fails  we  can  still  obtain 
the  Reynolds  stress  qw  (qw  = uw  + but  in  near  neutral  conditions 
w(w2  + v2  + u2)  ^ u3  so  wv2  < u3  . The  extra  term  < 1/50  u*  so  qw  = uw 


within  1 or  2%).  Note  also  that  f)  = 0 + ^ 1.005  U and  for  most  purposes 


either  may  be  used. 


1 
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Fig.  A. 2.  Deviation  from  cosine  response;  solid  line  6'  = 1 . 103  - 0.0047  0 
dashed  lines  have  slopes  of  - 0.0038  and  - 0.0056. 
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For  angles  of  attack  greater  than  10-20°  the  deviation  from  the  cosine 
> response  must  be  taken  into  account.  Since  6,  the  error  angle  if  the  axis 

of  the  v(  propeller  is  not  horizontal,  will  be  small,  the  formula  for  v, 
should  be  sufficiently  accurate.  The  observed  value  of  v',  v",  is  found 
to  be  smaller  than  expected  from  cosine  behaviour;  v"  = B'v.', . Figure  A. 2 
shows  values  from  four  factory  curves  (R.M.  Young  Co.)  as  well  as  values 
from  Hicks  (1972).  Hicks'  values  for  1.2  ms-1  are  somewhat  low  as  might 
be  expected  since  below  1 ms-1  the  output  is  lower  than  predicted  from 
the  linear  relation  between  output  voltage  and  wind  speed  for  higher 
speed  operation.  His  values  for  8 ms-1  are  somewhat  high,  exceeding  1.0 
for  o = 10-30°,  perhaps  because  of  wake  effects  on  the  vertical  shaft 
extension  which  was  added  to  make  the  geometry  more  symmetric  near  90°. 

In  the  range  of  40-70°  the  slope  is  similar  to  the  factory  curves.  We 
do  not  have  such  a shaft  extension  since  we  work  near  o = 60°.  For  any 
one  set  of  points  a linear  relation  in  the  range  35-70°  or  so  fits 
reasonably  well  with  some  variation  in  the  slope  and  intercept.  Gill 
(1975)  and  the  R.M.  Young  Co.  (personal  communication)  do  not  find  any  varia- 
tion of  e'  with  flow  speed  up  to  15  m/s  provided  it  is  sufficient  to  be  in  the 
linear  operating  range.  From  this  figure  we  take  p1  = 1.103  - 0.0047  o 
where  o is  the  angle  of  attack  in  degrees  and  this  formula  is  valid  for 
0 from  about  35-75°,  For  the  v,  propeller,  0 = a + 6 - tan-1(w/Q). 

Because  w and  Q are  negatively  correlated,  o is  more  likely  to  he 

< a + 5 than  > a + 6.  For  example,  if  w = -3o  =-0.18  Q (or  U the  mean 

w 

of  the  downstream  velocity  component  which  differs  from  Q by  about  l2«) 
and  taking  q = 0 (although  q is  likely  to  be  positive  reducing  the  angle 
difference)  0 = « + 5 + 10°.  For  w = 3a  and  q = -3o„  (where  = 1/10Q), 
o = a + 6 -14°.  Now  w and  q (or  u the  downstream  fluctuation)  have 
nearly  Gaussian  probability  distributions  so  the  chance  of  tan-1  w/Q  going 

outside  the  range  +10°  to  -14°  is  very  small  (less  than  1 in  500),  The 

tilt  angle,  5,  which  one  of  course  tries  to  make  as  small  as  possible, 
should  be  a few  degrees  at  worst.  In  practice  we  have  found  indicated  iS's 
to  oe  +5  - 6°  at  most.  Because  of  the  uncertainty  of  the  points  in  Figure  A. 2 
the  slope  which  is  taken  to  be  -0.0047  could  range  from  -0.0038  to  -0.0056 

as  shown  by  the  dashed  lines.  The  effect  of  these  possible  variations 

will  be  examined. 
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In  order  to  do  numerical  calculations  to  show  what  results  can  be  obtained 
we  shall  take  a = 60°.  In  practice,  of  course,  we  use  the  measured  value 
for  a particular  instrument  which  we  can  obtain  we  believe  to  rather  better 
than  1°. 

Converting  now  to  6 and  tan_1(w/Q)  in  radians  and  taking  a = 60°  for 
illustrative  purposes 

6'  = 0.821  [ 1 - 0.328  (6-tan-1  (w/Q))  ] . 

It  is  convenient  to  multiply  both  sides  of  the  equation  for  v"  by  1/0.821, 
that  is,  to  include  the  deviation  from  the  cosine  response  for  the  case 
w = 0 and  6=0  in  the  calibration,  then 

vn  = v2/0.821  = B [ Q cos(n/3  + 6)  + w sin(n/3  + 6)  ] 

where  & = 1 - 0.323  (6  - tan- 1 (w/Q) ) . In  the  case  a / 60°  we  use 

= B 1 (a)  ( 1 - (0.269/B *(«))( 6 - tan' 1 (w/Q)) ] 

and  take  v2  = v"/t*'  (a) . 

Now  w/Q  is  fairly  small  so  we  can  make  the  approximation  tan-1 (w/Q)  = w/Q. 

This  approximation  neglects  a term  of  order  (w/Q);1  but  here  we  shall  keep 
only  terms  to  second  order  in  small  quantities  (w/Q  and  6)  and  show  that 
they  give  sufficiently  accurate  results. 

The  tilt  angle  5 may  vary  slightly  because  of  cross-stream  fluctuations  v 
(from  existing  data  = U/10).  For  fluctuations  in  wind  direction,  4), 
p = tan-1(v/U)  = v/SJ.  6 will  have  its  largest  fluctuations  when  *p  is  such 
that  6 is  fluctuating  about  zero.  Taking  the  maximum  tilt  to  be  10°  (although  5° 
is  about  as  large  as  we  have  observed)  6 = 10°  s i ng.  10°  x p ~ 10°  v/U  and 
(T^)  2 = ° 10°  ov/U  « 1°.  Thus  in  the  worst  case  since  v and  hence  6 

will  be  approximately  Gaussian  we  need  to  consider  6 variations  up  to  *2° 
about  the  mean.  Ue  shall  look  at  the  effect  of  6 variations  after  considering 
the  case  when  6 is  constant. 

Now  we  have 

Vj  = Q cos6  + w sin6 

v , = [Q  cos  ( n/3  + 6 ) + w sin(  n/3  + 6)  ] • (1  - 0.328  6 + 0.328  w/Q} 


Taking  averages  in  time  for  a period  of  order  *4  hour  or  more,  w <<  Q 
so  we  assume  w = 0.  Then 

v = Q cos6 

v2  = Q [ cos  ( 11/3  + 6 )[  1 - 0.323  6 ] + 0.328  sin  ( n/3  + 6 ) vP/Q2  ] 

correct  to  second  order.  The  term  in  w2/Q2  is,  taking  o = 0.06  (J  and 
6=0,  0.001  while  the  other  term  is  0.500.  This  term  can  probably  be 
ignored  compared  to  other  errors  (in  calibration,  D.C.  drift,  etc.). 

The  procedure  we  use  is  to  remove  the  0 part  of  v2  by  taking 
v2-(v2/71)v1  and  then  solve  for  w.  When  we  ignore  the  w2/§2  term  we 
put  about  -0.001  Q and  hence  -0.001  u in  v,  - ( v.,/v,  )Vj . If  we  correlate 
this  with  u then  we  add  about  -0.001  u2  to  uw  or  about  0.006  uw 
for  an  error  of  about  ’a%.  For  our  four  blocks  of  about  13  minutes  each 
the  ratio  V^/v^  for  an  individual  block  varies  from  the  mean  for  the 
four  blocks  by  about  + 0.002  so  the  assumption  w = 0 should  introduce  a 
random  error  in  uw  of  less  than  IT.  Note  also  that  ignoring  this  term 
and  higher  order  terms  may  give  a w with  a non-zero  mean  so  it  is 
important  to  remove  the  mean  in  w and/or  in  u or  integrate  the 
cospectrum  leaving  out  the  zero  frequency  value.  If  we  calculate  uw  as 
(U  + u)  x (w  - .001  (U  + u) ) we  qet  uw  0.001  0 -0.001  7 or  about 
uw  (1  + 0.6)  which  is  a very  seriou  error.  Thus  here,  as  in  any  Reynolds 
stress  measurement,  it  is  essential  to  remove  the  means.  Actually  the 
calculation  is  a little  more  complicated  because  with  v./V,  = 0.501  the 
calculated  6 is  not  zero  but  -0.001  radians.  Then  u is  calculated  from 
v,  - Q - w tan6  and  is  actual  ly  u + . 001  w and  (w  - .. 001  u)  • [u  + . 001  w)  = uw  - . 001 
(u^-w7)  = uw  ( 1 + .004)  so  the  error  is  about  \ *. 

Clearly  for  this  procedure  of  removing  the  0 part  of  v,  to  work  one 
must  have  good  D.C.  stability.  We  have  been  very  careful  about  D.C. 
stability  in  designing  our  circuits  using  low  temperature  coefficient 
zener  diodes  as  references  for  our  power  supplies  and  offset  voltages. 

Before  discussing  the  effects  of  offset  errors  further,  because  they 
produce  apparent  6's  which  must  be  taken  into  account,  let  us  first  show 
examples  when  there  are  no  errors. 

Neglecting  the  w-Vq2  term  for  the  moment  we  have 
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v^/Vj  = cos(n/3  + 6)  (1  - 0.328  6 )/cos6 

= (1  - 0.328  5 ) (cos  n/3  cos5  - sin  ll/3  sin5)/cos6 

= (1  - 0.328  5 ) (1  -/3  tan5) 

Since  5 is  small,  take  tan  6=5  which  will  give  an  expression  correct  to 
order  63.  Then  2v,/Vj  = [1  - 2.060  6 + 0.568  62]  which  allows  us  to  calcu- 
late 5 from  the  observations  of  v’j  and  v2.  The  quadratic  formula  gives 
6 = 1.813  ± ✓1.528  + 3.521(72/1^).  Since  6 is  clearly  less  than  1.813  radians 
we  see  that  it  is  the  root  involving  the  -/”  that  is  required. 

Now  v,  = Qcos(n/3  + 6)  (1  - 0.328  6) 

+ w[  sin(n/3  + 6)  (1  - 0.328  6)  + 0. 328  cos(n/3  + 6)  ] 

• I 

+ 0.328  w2  si n (n/3  + 6 ) /Q 

v(  = Q cos6  + w sin6  > 

so  v , - (v2/7j  )Vj  = w [ [ sin(n/3  + 6)  - cos ( n/3  + 6)  tan6]  [1  - 0.328  6 ] 

+ 0.328  cos( n/3  + 6 ) ] + 0.328  w2  sin(n/3  + 6)/Q 
- 0.328  sin(n/3  + 5)  (w2/Q2)  (Q  cos6+w  sin6)/cos6  . 

Using  sin(n/3  + 6)  cos5  - cos(n/3  + 6)  sin6  = sin  n/3  we  have 

v2  - ( v^/Vj  )vj  = w[ (sin(n/3yco$6) ( 1 - 0.328  6 ) + 0.328  cos(n/3  + 5)] 

+ 0.328  w2  sin(n/3  + 6)/Q 

- 0.328  sin(n/3  + 6)(w2/Q2)(Q  cos6  + w sin6)/cos6. 

Now  the  cos(n/3  + 6)  w term  arises  from  the  ? Q cos(n/3  + 6)  term.  Thus  even 

: : 

to  first  order  it  is  not  sufficient  to  simply  correct  for  the  deviation 
from  cosine  response  for  the  w = 0 angle  of  attack.  Using  the  linear 
approximation  and  neglecting  this  term  gives  w's  almost  20%  too  large. 

The  uncertainty  in  the  coefficient  0.328  will  lead  to  some  error  hut  it 

is  only  about  3%  as  we  shall  show  later.  The  final  term  can  be  neglected 

without  serious  error  provided  we  remove  the  mean  from  u and/or  w before 

calculating  uw  as  noted  before.  The  term  in  w2  is  kept  because  w2  ♦ 

may  be  several  times  w2.  However,  with  negligible  error  we  can  replace  Q 

with  Vj/cos6  in  this  term.  So  the  approximation  we  shall  use  is 

j 
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v2-(v2/v1)v1  = B w + A w2/Vj 

where  Vj/Vj  observed  is  cos(n/3  + fi)(l  - 0.328  6 )/cos6  + 0.328 
• sin(n/3+5)  w2/Q2 

B = sin(n/3)(l  - 0.328  6 )/cos6  + 0.328  cos(n/3  + 6) 

A = 0.328  sin(n/3+6)  coss 

and  6 is  calculated  from  1.813  - /1.528  + 3.521  (v^/v^) 

First,  let  us  consider  some  examples  where  we  calculate  v:  and  v2  without 
any  approximations  using  assumed  values  of  Q,  w,  and  6.  Then  we  calcu- 
late w from  the  approximate  formula,  is  the  solution  of 

v2  - (v2/v'1)v1  = B w^  ; wq  is  the  solution  of  the  quadratic  where  with 
B>0  the  correct  root  involves  +/b2  - 4ac  in  the  equation  aw2  + bw  + c = 0. 

We  take  Q = 10  ms-1  o = 0. 1 $ = 1 ms-1  and  a = 0. 06  Q = 0. 6 ms-1 . Since 

H w 

u (and  hence  Q)  are  negatively  correlated  with  w we  take  ±2 ow  with  + 2oq 
and  ± 3a  with  +3 a . These  include  about  96%  and  99.7%  of  cases  for 

w q 

Gaussian  variables  for  CD  equal  to  1.6  xl0‘3.  When  CD  = 2.1  xlO"3,  then 

the  ± 3a  case  given  includes  99%  of  cases  provided  a /u*  and  o /u+  remain 

w g 

the  same  (as  expected  since  observations  over  land  give  similar  values). 


Take  6 = 

10°: 

v2/ v j observed  = 

0.3285;  6 

calculated 

is  10.02°. 

Calculated 

A and  B 

are 

0.3036  and 

0.9411  (for  6 = 10 

°A  & B are 

0.3035  and 

0.9412). 

w 

Q 

vi 

V2 

w* 

w„  w„ 
q c 

, .OOllv, 

1 B 

1.8 

7 

7.206 

4.189 

1.936 

1.792 

1.800 

1.2 

8 

8.087 

3.831 

1.248 

1.191 

1.201 

0 

10 

9.848 

3.224 

-0.012 

-0.012 

0.000 

-1.2 

12 

11.609 

2.709 

-1.174 

-1.215  - 

1.201 

-1.8 

13 

12.490 

2.473 

-1.733 

-1.818  - 

1.803 
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w,  = w+fA/BVjJw2  and  shows  the  expected  distortion.  Later  we  shall  show 
that  the  effect  of  ignoring  the  quadratic  term  is  reasonably  small,  w 


differs  from  w mainly  by  - ( 0 . 001 1/B) v 1 which  comes  from  neglecting  the 
w2/Q2  term,  q or  u calculated  from  v^os 6 - w tanfi  will  be  in  error  by 

10-4 

in 


u-4xl0‘4w  which  is  negligible.  As  noted  before,  the  10-3  u term 
w calculated  gives  an  error  of  about  + '*%  in  |uw|. 


Take  6 = 0 : v ,/v1 


observed  = 0.5010; 


6 calculated  is  - 0.06° 

for  <5  =0° 


B = 1.031 
B = 1.030 


A = 0.2839 
A = 0.2841 


w 

Q 

vi 

V2 

wq  w 

+ - 
q 

.0010  v, 

B 

1.8 

7 

7. 

5.476  1 

.911 

1 

.785 

l 

.792 

1.2 

8 

8 

5.285  1 

.239 

1 

.191 

l 

.199 

0 

10 

10 

5.000  -0 

.010 

-0 

.010 

0 

.000 

-1.2 

12 

12. 

4.799  -1 

.177 

-1 

.211 

-l 

.199 

-1.3 

13 

13 

4.718  -1 

.741 

-1 

.811 

-l 

.798 

error 

is  - 

0.001 v. 

. The  q 

error  • 

is 

lO^u  + 

10' 

‘3  w and  uw 

-3(ir'  - w 

2)  < 

*a  %• 

» 5 = - 

10°: 

1 

observed  = i 

3.6910; 

<5  calculated  is  - 

10. 13° 

B = 

1 

.142 

A = 0.2469 

for  6 = -10° 

B = 

1 

.141 

A = 0.2474 

w 

Q 

V. 

V 

w 

w 

w 

. 0009  v \ 

...  ■ - ■ — » 

l 

l 

q 

1 

q b 

1.8 

7 

6.581 

6.700 

1.835 

1.781 

1.786 

1.2 

8 

7.670 

6.705 

1.230 

1.190 

1.196 

0 

10 

9.848 

6.796 

-0.008 

-0.008 

0.000 

-1.2 

12 

12.026 

6.961 

-1.181 

-1.207 

-1.198 

-1.8 

13 

13.115 

7.062 

-1.752 

-1.805 

-1.795 

The  q error  is  - 5 x 10'4  u + 2 x 10-3  w.  Error  in  uw  calculated  v'*  %. 
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Clearly  given  a steady  value  of  6,  an  accurate  a,  good  D.C.  stability, 
accurate  calibrations  and  an  accurate  formula  for  t<’,  we  can  calculate 
w with  good  accuracy  and  with  negligible  contamination  with  u,  the 
downstream  velocity  component.  Removal  of  u in  the  calculated  w is 
essential  to  avoid  large  errors  in  uw  in  any  method  of  obtaining  uw 
and  the  procedure  is  designed  specifically  to  do  this  removal  as  well  as 
possible.  Note  that  even  at  5 ms-1  and  6 = 10° , v,  > 1 ms'1  so  that  we 
stay  in  the  linear  range  of  the  propeller's  output,  the  original  intent  of 
using  a = 60°  rather  than  90°.  At  higher  turbulence  levels  (e.g.  over 
fairly  rough  land)  where  6 might  become  as  large  as  10°,  it  would  probably 
be  best  to  use  a = 50°  . A smaller  a would  decrease  the  relative 
importance  of  the  guadratic  term  somewhat  but  increase  the  uncertainty 
in  B and  6 because  of  the  uncertainty  of  the  slope  in 


Let  us  now  look  at  how  various  errors  affect  the  calculated  w.  Since  w 
is  also  used  in  calculating  u (or  g)  from  vl  these  errors  may  propagate 
into  u and  the  calculated  uw  so  the  effect  on  uw  must  be  examined 
as  well . w is  obtained  from  a w2  + b w + c = 0 where  c = - (v2  - (vV^i  )vi ) ♦ 
b = B and  a = A/Vj.  Now  |4ac|  <<  b so  we  can  expand  w = (-b  + /b^"-  4ac)/?a. 
to  give  w = - c/b  - (a/b)(c/b)2  in  order  to  do  our  error  analysis.  Note 
that  this  expansion  is  equivalent  to  using  w ^ in  the  quadratic  term  to 
get  an  approximate  solution  for  w 


w^,  = c/b  and  w = w , - (a/b)w-'  . 

Calculate  w’  = w^,  - (a/b)w^,2  = w + (a/b)w2 -(a/bXw2  + (2a/b) w3  + ( a2/b2)w4). 
w 1 = w + 0(w3)  so  w’  is  also  a second  order  accurate  solution  and  could  be 
used.  In  practice  w^  the  solution  to  the  quadratic  is  better  containing 
smaller  w3  terms  than  w’.  However,  in  analyzing  for  the  error,  sw, 
due  to  variations  in  a,  b,  c the  second  order  accurate  w’  solution  is 


easier  to  use  and  should  be  of  sufficient  accuracy. 
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The  first  two  terms  came  from  shifts  in  w due  to  c and  b changes  while 

< 

the  final  three  are  associated  with  changes  to  the  'correction1  term  -(a/h)w,2. 
As  this  term  is  quite  small  ( < 10%  of  w.)  the  three  final  terms  will  generally 
be  of  secondary  importance.  If  we  calculate  Jw  from  (?w/ aa)sa  + (a  w/?  b^b  + 
(aw/acjiic  then  the  terms  above  are  the  largest  terms  in  an  expansion  of 


this  more  exact  expression. 


Errors  in  a 

Suppose  a is  actually  62°  instead  of  60°.  v2  will  then  be  smaller  than 
expected  because  we  have  included  b'(60)  in  the  calibration,  not  3 1 ( 62 ) ; 
the  factor  is  0.9886.  With  6 = 0,  v2/7,  observed  is  0.9886  cos(62)  + 

0.328  sin(62)  w^/Q2  = 0.4651.  6 calculated  is  1.96°  because  in  the  formula 

used  to  calculate  6 the  1 - 0.328  6 factor  accounts  almost  exactly  for  the 
6,(62Vb'(60)  factor  (if  the  w^/Q2  term  is  subtracted  from  72/Vj  then  6 
calculated  is  2.02°).  Thus  any  error  in  a appears  as  an  apparent  tilt 
angle  6. 

A should  be  0.9886-  (.328/0.9886)  sin62°=  0.2896, 

calculated  to  be  0.328  sin(61.96°)  cos(1.96°)  = 0 . 2893; 

B should  be  0.9886  sin62°+ 0.328  cos62°=  1.027, 
calculated  to  be  1.011. 

There  is  no  error  in  c as  it  is  taken  as  given;  there  is  a small  error 
in  a (6a/a  = 0.015)  but  as  it  is  in  a second  order  correction  the  error  in 
w is  2x10" 3 at  most  and  can  be  ignored.  The  main  error  is  in  b ( =B). 
Since  B used  is  too  small,  w will  be  too  large  by  about  1?$. 


w 

Q 

vi 

V2 

w 

q 

V-Y 

1.8 

7 

7 

5.222 

1.945 

1.311 

1.818 

1.2 

8 

8 

4.995 

1.261 

1.208 

1.216 

0 

10 

10 

4.641 

-0.010 

-0.010 

0.000 

-1.2 

12 

12 

4.372 

-1.196 

-1.232 

-1.220 

-1.8 

13 

13 

4.258 

-1.763 

-1.843 

-1.830 

w is  a little  large  as  predicted  and  contains  -0.001  u.  However,  because 

we  take  u to  be  the  fluctuating  part  of  Vj/cosfi  - w tan6  and  v]  actually 
has  no  w in  it  (the  true  6 = 0)  we  get  1.0006  u - 0.034  w.  uw  calculated 
(1.015w- 0.001  u)(1.0006u  - 0.034  w)  = 1.016  uw  - .001  u7  - 0.035  w7.  With 

0.,  = 1.5  u*,  o = 2.5  u*,  | uw]  calculated  is  10%  too  large.  Thus  while  we 

w u 

have  gotten  the  u out  of  w pretty  well  we  have  left  some  w in  the 
calculated  u.  This  problem  is  not  peculiar  to  this  system.  An  error  in 
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the  angle  between  the  axes  along  which  components  are  resolved  will  lead 
to  the  same  problem  in  any  system  and  the  error  will  be  similar  assuming 
one  is  careful  to  use  w = 0 to  remove  any  u signal  from  the  calculated 
w.  The  error  is  a little  worse  here  because  the  w calibration  is  also 
affected  a little  (^1.5%).  For  a one  degree  error  in  a the  uw  error  is 
about  5%;  for  a system  for  which  only  the  w in  u causes  an  error  (e.g. 
a sonic  anemometer)  the  error  would  be  4*.  We  believe  we  can  measure  a 
to  rather  better  than  1°  but  to  be  on  the  safe  side  we  must  consider  there 
to  be  a possible  error  of  up  to  ± 5%  from  a errors.  Since  0 is  not  appre- 
ciably affected  the  error  in  Cp  is  also  similar. 

Offset  errors 

We  have  been  careful  to  make  these  very  small  as  noted  earlier,  and  we  expect 
them  to  be  similar  in  both  v:  and  v,  in  % terms  so  that  they  cancel.  At 
5 m/s  an  offset  error  corresponding  to  ± 2 bits  (i.e.  ± 5 mv)  is  about 
1/200  in  Vj  and  1/150  in  v2.  Suppose  the  worst,  that  these  offsets  occur 
with  opposite  signs.  With  Q = 5 ms-1  and  6=0  for  Vj  we  add  0.025  ms-1 
to  the  values,  and  for  v2  we  add  -0.015  ms-1.  v2/vj  = 2.485/5.025 + .001  = .4955 
instead  of  0.5010,  6 calculated  = 0.25°.  As  well  as  a small  offset  error  in 
Wq,  which  we  ignore  since  means  are  removed,  we  get  w + 0.0045  u.  A cal- 
culated is  0.2848,  correct  A is  0.2841.  D calculated  is  1.028;  B cor- 
rect is  1.030.  The  main  error  is  due  to  adding  some  u to  w.  uw  cal- 
culated is  (w  + 0.0045  u)  (u  - .0044  w)  = uw  + .0045  (u2  - w2) , so  |uw|  is 
about  2%  too  small.  U2  is  about  1%  high  and  Cp  about  3%  too  small.  If  the 
errors  are  opposite  v.,/^  = 0.5065,  w calculated  = w-.065  u,  6 = -0.36°, 
u calculated  = u+.0063  w.  Error  is  -0. 0065{u2  - w2) ; fuv7|  is  about  2.5% 
too  large  and  Cp  about  3.5%  too  large.  Thus  at  5 ms-1  the  error  is  about  2.5%, 
at  10  ms-1  about  1%  in  uw  and  perhaps  2%  in  Cg,  and  lower  at  higher  speeds. 
Note  that  the  apparent  6 is  such  as  to  add  a little  w to  the  calculated 
u which  helps  to  reduce  the  error.  In  fact  the  offset  errors  due  to  thermal 
drift  are  in  the  same  direction  so  the  actual  errors  are  \ or  less  than 
these  worst  case  values. 


Calibration  errors 

For  one  propeller  relative  to  the  other,  the  calibrations  are  probably  accurate 
to  about  ± 1%.  Overall*  based  on  our  own  checks  amongst  others,  they  should 
certainly  be  better  than  ± 2%.  If  both  propellers  are  2%  high  or  low  then 
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the  effect  cancels  in  calculating  5,  A,  B;  w and  u will  both  be  2?  high 
or  low  together;  uw”  will  be  ± 4?  but  in  calculating  Cp  the  error  cancels 
out  and  Cp  is  unaffected.  If  v,  is  2%  high  while  v.,  is  2 % low  and 
5 true  = 0,  then  v,/Vj  = 0.4814.  In  taking  v -(v./VjJVj  we  are  still 
left  with  - 0.001  u;  using  the  observed  ratio  still  gets  the  u out  of 
w in  the  same  way  as  before,  w is  of  course  2%  too  small  due  to  the  cali- 
bration error,  but  B calculated  is  1%  low  so  w is  only  1%  low,  u is 
2%  too  big,  6 apparent  is  1.04°,  so  u apparent  is  1.02u-.018w  and 
uw  as  calculated  is  1.01  uw  - . 001  u2  - . 018  w2  and  is  about  5? 
too  large  in  size;  again  the  effects  on  Cp  tend  to  cancel  and  it  is  only 
1?  too  large. 

If  is  2*  low  while  v.  is  2$  high  and  6 true  = 0,  then  v./Vj  = 0.5214, 

5 calculated  is  -1.18°;B  is  1%  high  and  w is  1%  high;  u is  2%  low, 
but  has  0.021  w in  it  so  uw  calculated  is  ( . 98  u + .021  w) ( 1 .01  w - .001  u) 

= 0.99  uw  + 0.021  w2  - .001  u' . uw  is  about  5"  too  small  in  size  but  since 
U is  2%  low,  Cq  is  only  about  1"  low. 

Errors  in  the  6'(o)  relation 

From  Fig.  A. 2 it  is  clear  that  some  variation  in  this  function  is  possible. 
To  examine  the  effects  let  us  assume  that  the  relation  used  so  far  is 
correct  and  see  what  happens  when  we  use  some  other  B‘(e)  in  calculating 
w,  uw  and  Cp.  Suppose  we  think  that  ;; ' (60)  the  intercept  is  0.84  instead 
of  0.821,  then  the  calculated  v ‘s  will  be  too  small  by  0.821/0.84  = 0.977 
and  so  will  v,/7j  by  the  same  amount.  For  true  s's  of  10,  0,  -10  the 
calculated  iS's  are:  10.49,  0.59,  -9.32.  v,-(v,/v,)v[  will  be  2.3?  low 

and  will  contain  - 0.001  u as  before.  Because  of  the  apparent  5,  B will 
also  be  too  low  by  about  0.6*  and  A will  be  about  0.6"  high.  The  change 
in  A can  be  ignored  since  it  occurs  in  a small  correction  term.  Thus 
| w | calculated  will  be  about  1.7*  too  small.  At  the  same  time,  we  intro- 
duce some  w into  the  calculated  u (-.9,  -1.0  and  -1.2*  for  6 ~ 10, 

0,  -10°).  So  | uw | calculated  has  an  error  of  -1.7  from  small  w,  +0.6  5. 
from  u in  w and  +1.9  to  2.5:  from  w in  u due  to  the  apparent  <S. 

Thus  ) uw]  is  0.8  to  1.5?  too  large  and  so  is  C,,.  If  the  assumed  inter- 
cept is  0.80  then  |uwl  and  Cp  are  almost  unchanged  because  the  effect 
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of  the  - 0.001 u in  w cancels  the  intercept  change.  Our  original 
calculations  were  based  on  the  curve  with  the  0.82  intercept.  When  the 
other  data  are  plotted  it  appears  that  0.81  might  be  better  in  which  case 
our  | uw | ' s are  about  \%  too  small,  neglecting  the  - 0.001  u in  the 
calculated  w.  However,  neglecting  the  w2/Q2  term  in  v,/Vj  tends 
to  make  |uw|  about  '2%  too  large  and  compensates.  The  effect  of  this 
possible  error  in  the  intercept  is  about  + 1%  for  intercept  variations 
of  ± 0.04  and  | 5 | < 5°. 

With  the  intercept  unchanged,  variations  in  the  slope  used  in  the  calcu- 
lations do  not  change  Vo/vj  but  do  channe  A,  8 and  the  apparent  &. 

With  the  slope  variations  shown  in  Figure  A. 2 the  0.328  coefficient  in 
the  equations  may  vary  from  0.265  to  0.391.  A then  also  varies  by 
± 20%.  The  effect  of  using  the  wrong  A is  then  to  leave  about  20% 
of  (A/Bv^w2  or  second  harmonic  distortion  in  the  calculated  w.  As 
we  shall  show  later,  the  effect  of  ignoring  the  whole  second  harmonic 
(w2)  term  is  small  so  if  20%  of  it  is  left  in  it  should  not  matter. 

Suppose  the  correct  coefficient  is  0.328  but  we  assume  it  to  be  0.265, 
then  the  5 equation  becomes  1 - (/J+  0.265)5  + 0. 265/3  S2  = 1 - 1.997  5 + 

0.459  52.  For  an  actual  5 = 10°,  5 calculated  is  10.38°,  B calculated 
is  1.5%  low,  making  jwj  calculated  1.5%  too  large,  u calculated  is 
actually  u-  .007  w so  |uw|  calculated  is  3%  too  large  and  so  is 
(actually  3.5%  considering  the  u in  w effect).  At  5 = 0 only  the 
effect  on  B needs  to  be  considered;  it  is  3%  low  and  |w|  and  |uw| 
are  3%  high.  For  6 = -10°,  5 calculated  is  -10.63°,  B is  4%  low  and 
w,  4%  high,  but  u calculated  is  u + O.Ollw  so  |Iiwj  is  in  error  by 
+ 2%. 

If  we  assume  the  slope  coefficient  is  0.391  instead  of  0.328,  the  5 
equation  is  1-2.1235+  0.677  52.  For  5 = 10°,  5 calculated  is  9.79°, 

B is  1.5%  high  and  |w|  1.5%  low;  u calculated  is  u + 0.004  w,  so 
|uw|  is  2.5%  too  snail  from  the  slope  effect  and  the  - .001  u in  w 
reduces  this  to  < 2%.  At  5 = 0,  B is  3%  high  and  |w|  3%  low;  |uw| 

is  2.5%  low  including  the  effect  of  the  - 0.001  u in  w.  For  5=  -10°  #j 

5 calculated  is  -9.78,  R is  4%  high  and  |w|  4°  low;  u calculated  is 

u-  0.004  w so  | uw ' is  2.5%  low  including  the  u in  w effect. 

, 1 
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If  the  slope  and  intercept  both  change  the  effects  should  be  additive. 

From  Fig.  A. 2 a low  slope  would  tend  to  go  with  a high  intercept  and 
vice  versa.  At  s = 0,  we  get  a total  error  of + ^%  from  the  intercept, 

+3%  from  slope,  + *s%  from  the  u in  w for  a total  of  4%.  By  calcula- 
tion, | w | is  too  small  by  2.3%  from  the  change  in  v,,  -(72/71)v1,  6 

calculated  is  0.61°,  B is  3.5%  too  small,  so  w is  1.2%  too  large  and 
contains  - .001u.  u calculated  is  u-.Ollw  so  |uw|  is  +1.2%+  0.6% 

+ 2.2  = 4%  too  large,  which  agrees.  For  a high  slope  and  low  intercept  at 
6 = 0,  |uw|  is  - H%  from  the  intercept,  + %%  from  the  u in  w and 
-3%  from  the  slope,  for  a total  of  -3%. 

Thus,  overall  the  error  in  e'(e)  should  lead  to  an  error  of  about  ± 3%. 

As  the  slope  for  any  given  data  set  (where  the  experimental  error  does 
not  cause  a lot  of  scatter,  i.e.  beyond  about  70°)  seems  to  be  well 
within  the  limits  which  we  have  used  to  calculate  errors,  the  error  is 
probably  less  than  + 3%. 

Errors  due  to  fluctuations  in  6 

6 = 6 cos U-<j>  ) where  * is  the  wind  direction  and  6 has  its  maximum  6 

O O 0 

when  <j>  = 4>q.  Taking  <j>  to  be  in  the  direction  of  U the  mean  downstream 
component 

4>  -7  = sin"1  v/Q;  sin(<f>-7)  = v/Q 

COSU-*)  = 1 - v2/2Q2  + 0(v4/Q4)  and  the  term  0(v4/Q4)  <^0.001 

99.7%  of  the  time  if  v is  Gaussian,  so  we  neglect  it. 

6 = 6q{cos(4>  - ?)  cosUq -7)  + sin(4>-7)  sin(<j>p-7)) 

= 6o{cosUo  -?)  ( 1 - v2/2Q2)  + sinUo -♦)(v/Q)} 

7 = 6o(cosUo-7)  (l-77/2Q2  + higher  order  terms) 

Fluctuations  in  6 = 6 - "6  = 6 {cos($  - <t>)(VTT/2Q2  - v2/2Q2)  + sin(<}>  -7)(v/Q)l. 

0 0 0 _ 

When  4>  -$  = 0,  7=6  (1  - .005)  so  the  difference  between  6 and  6 cos(<j>  -4') 
0 0 0 0 

is  about  0.05°  at  worst. 
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Note  that  the  6 fluctuations  have  both  a linear  and  a quadratic  term  in  v/Q. 


Consider  the  Q part  of  v,  with  6 = 7 + 6' 

Q(cos  7cos  <5 1 - sin  6" sin  6 ' ) 

0 cos  61  = Q [ 1 - 62  sin2 ^ - 4>)  v2/202]  to  second  order 

= Q cos  o 


Q sinks'  = 0 to  second  order 


We  see  that  although  Q and  <s ' are  correlated  it  is  at  such  high  order 
that  we  may  ignore  it. 


vo  = Q [ cos(n/3  + 6)  cos  6 1 - sin(n/3  + s i n 6 ' ] [1  - 0.328(<s  + 6 ' ) 1 


+ f(w) 


v = 0 cos  (n/3  + 6)  cos  a,  (1  - 0.328  6) 


- Q sin  (n/3 + 6)  • 0.328otf2  + 0 (vr/0) 


but  the  second  term  is  0.0001  Q which  is  negligible  so  we  get  the  correct 
v2/Vj  and  6 even  if  6 is  fluctuating.  In  solving  for  w,  7 /v  , 6, 

A and  B are  the  same  as  if  6 were  fixed  but  as  6 chanqes,  v„  - (v„/v.)v, 
will  contain  extra  fluctuations.  Consider  first  the  part  involving  Q 


v = Q(cos  6 cos  6' - sin  <s  sin  5 ' ) 


v , = Q{cos(n/3  + 6)  cos  6 1 - s i n ( n/3  + 6)  sin  tf}{  1 - 0.328(6  + 6' ) } 


v?  - (Vj/VjJv,  = Q{cos(n/3  + J)  cos  6 ' • (-0.328  6') 

- si n ( 11/3  +7)  sin  6'  (1  - 0.328(7+6')) 
+ v_,/v1  sin  6 sin  6 1 - 0.001) 


Now  this  is  to  be  correlated  with  u or  q but  they  are  unrelated  to  6' 
to  second  order.  The  largest  term  (other  than  the  - 0.001  u?  associated 
with  neglecting  the  vT/Q  in  v.)  which  will  be  added  to  uw  is 
0. 328  sin  (n/3  + 6)  6'2.  u?  lO'NT'*  which  is  completely  negligible.  The 
fluctuations  in  w are  of  order  ± 0.04  Q for  6 = 10°  and  a.  = 1° 

O 6 

which  are  not  so  small  but  the  average  effect  is  very  small.  Likewise, 
there  will  be  a variation  in  the  calculated  w involving  w6*  and  w 6* 
terms  but  these  will  also  have  a completely  negligible  effect  on  the 
calculated  uw. 
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Non-cosine  response  of  the  V]  propeller 

Up  to  about  10°  the  deviation  from  cosine  response  is  undetectable  within 
experimental  error.  With  |<$|  < 5°  the  angle  of  attack  stays  within  10° 
almost  all  the  time.  However,  with  larger  tilt  errors  of  10-  15°  the 
effect  may  be  of  importance.  From  Fig.  A. 2 near  e = 15°,  g'  = 1 - 0.038  0 
with  e in  radians,  is  a possible  fit  and  v2  will  be 

Vj  = (Q  cos  6 + w sin  6)  (1  - 0.038  | 5 - w/Q| ) 

Suppose  this  relation  holds  even  when  6=0  and  consider  the  situation 
when  6=0 

Vj  = 0 - 0.038  | w| 

7 ' = Q - 0.038  JwJ 

Now  w is  approximately  Gaussian,  so  JvTf  = a //rf  - 0.06  Q//H  and 

w 

Vj  = 0.999  ().  72/7x  is  now  0.5015  and  6 calculated  is  -0.08°.  The 
effect  on  A and  B is  about  0.1%  which  can  be  ignored,  w calculated 
is  w - 0.0015  u + 0.019  |w|  while  u calculated  is  u + 0.0015  w - 0.038  |w| 
and  uw  calculated  is  "uw  - .0015(7^ -"w5)  - .0007  w” 5 , since  | w | 2 = w2, 
w|w|  = 0 and  u|w|  is  small  and  assumed  zero.  The  error  in  |uw|  is 
about  0.8%  and  the  error  in  Cp  = + 1%.  If  | 5 | is  larger  than  10°  or  so 
then  |fi-w/Q|  is  6 - w/Q  for  6 > 0 and  - 6 + w/Q  for  6 < 0.  There  is 
rather  more  w in  Vj  than  included  in  the  formulae  given  earlier  and 
it  may  cause  fairly  large  errors. 

Consider  6 = 10°,  then  taking 

Vj  = (Q  cos  6 + w sin  6)  (1  - 0.038(6  - w/Q)) 

Vj  = Q cos  6 (1  - 0.038  6)  + 0.038  sin  6 w2/Q 

The  second  term  is  » 2xl0-5  Q and  may  be  ignored.  = Q cos  6 x 0.0.993 

and  is  0.7%  smaller  than  if  the  possible  non-cosine  effect  is  ignored. 
v2/Vj  is  larger  and  6 calculated  is  9.89°.  The  error  in  B is  about 
0.1%  and  in  A even  less.  v2-(v2/7J)vi  is  as  calculated  in  the  simpler 
formula  except  for  the  extra  w term  which  is  - 0.012 w and  the  calculated 
w is  1.2%  too  small,  u as  calculated  is  0.993  u+  0.04w  and  uw  calculated 
is  (0.938w-  0.001  u)(0.993u  + 0.04  w)  = 0.981  uw  - .001  u2  +0.04  w2,  or  about 
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10%  too  smal'  in  size;  CD  is  9%  too  small. 

When  6 =-10°,  1 6 - w/Q | = 10°+w/Q  Vj  is  0.993  Qcos6  as  for  6 = +10°. 

In  v2-(v2-71)v1  the  extra  term  is  +0.026w.  6 calculated  is  -10.38° 

making  B 0.4%  too  large,  so  the  net  effect  is  that-w  is  2.2%  too  large, 
u as  calculated  is  0.993  u - .030  w.  uw  as  calculated  is  (1.022w  - .001u) 
(0.993  u - .030w)  which  is  about  9%  too  large;  is  about  10%  too 
large. 

From  the  experimental  data  it  is  difficult  to  determine  how  important  this 
possible  error  may  be.  For  small  |s|,  say  less  than  2°,  not  only  is  the 
deviation  from  cosine  response  likely  to  be  smaller,  but  also  the  contri- 
bution from  the  extra  w term  in  tends  to  cancel  as  it  does  when 
6=0  where  it  is  a |w|  term.  If  one  has  a large  body  of  data  then  one 
can  look  for  systematic  differences  between  6 = 0,  6 = 4-5°  and-6  = 4-5° 
in  the  data  set  and  we  shall  attempt  to  do  so  in  the  analysis  of  the  data  we 
have  obtained.  Some  further  wind  tunnel  work  would  probably  be  helpful 
too  in  the  range  e = 0 to  ± 20°. 

The  importance  of  the  quadratic  term 

It  is  a fairly  simple  matter  to  solve  for  wq.  However,  if  the  linear 
approximation  is  adequate,  operations  could  be  done  on  the  Fourier  coeffi- 
cients of  Vj  and  v2  which  have  to  be  found  to  correct  for  instrument 
response  and  then  inverse  transformed  to  calculate  w . 

w^  = w + aw2,  where  a = A/(BQcos6) 

wfc2  = w2  + ^aw3  + a2wU 

Put  Q = U + u 

w ^ = w2+  [2A/(B  cos  6 U)][w3  - uw3/U  1 + (A/(B  cos  6 ) )2  w4/U2 

w3  will  be  small  compared  to  ow3  and  w4  -v  3 aj4.  A/B  cos  6 is  0.33 
to  0.22  for  6 = +10  to  -10°.  The  w3  term  * 2/3  (o/U)  • (o  2)  • (w^/o  3). 

^ — ^ w w w 

ow/lT  <v»  0.06  and  w3/o  3 will  certainly  be  less  than  1/10  .’.  the  w3  term 
< 0.004  w7  and  can  be  neglected.  Even  if  u and  w3  are  perfectly  correlated 


i 


55 


the  uw3  term  < 0.004w2  . The  w4  term  < 0.0004  w2.  Thus,  w?2  = w2 

with  less  than  1%  error.  The  spectra  of  the  two  should  also  be  similar 

except  perhaps  where  the  w spectrum  is  very  small. 

Consider  ue  = Vj/cos  6 - w£  tan  6 - (U  + (v2/2U)) 

= u - A/(B  cos  6)-tan  6*(w2/U)  [1-u/U]  + higher  order  terms 
u?2  = u2  - (2A  tan  6 )/ ( B cos  6)  (uw2/U)  + higher  order  terms 

For  | £ | < 10°  the  second  term  is  s 0.005  u^"  even  if  u and  w2  are 

perfectly  correlated.  Thus  is  a good  approximation  for  u. 

Finally,  consider  tTjw? 

= (w  +(\/(B  cos  6))  (w2/Q) ) (u -(J A tan  6)/(B  cos  <5  Q))w2) 

= uw  +[A/(Bcos  5)]  [uw2/Q  - fo3  tan  6)/Q] 

The  largest  error  term  is  A/(Bcos6  ‘(Duw^.  For  perfect  correlation  this 
term  would  be  about  7 % of  |uw|  but  the  correlation  will  be  much  less  than 
1 so  the  error  is  probably  negligible. 

At  present  we  are  calculating  wq  not  w ^ but  we  plan  to  compare  results 
from  calculations  using  both  with  real  data  to  check  that  the  use  of  w{ 
is  adequate  as  it  appears  to  be. 

Because  of  slope  errors  about  20%  of  the  (A/B  cos  <s) (w2/Q)  term  may  be  left 
in  even  when  we  calculate  w^  but  as  the  effect  of  the  whole  term  is  small 
leaving  in  20%  of  it  can  be  ignored. 
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